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ABSTRACT 



Aims. The dust emission spectrum and the brightness profile of passively heated condensed cores is analyzed in relation 
to their astrophysical environment. The model is used to study systematically the radiative transfer effects on essential 
parameters such as the dust emissivity, dust temperature, and luminosity of the cores and to derive uncertainties in 
typical estimates of the IR flux and size. 

Methods. The cores are modeled as critically stable self-gravitating spheres embedded at the center of self-gravitating 
filaments that are assumed to be either spherical or cylindrical in shape. The filaments are heated by an isotropic 
interstellar radiation field (ISRF). The calculations are based on a physical dust model of stochastically heated grains 
of diffuse interstellar dust. The spectral energy distribution (SED) of the cores is calculated using a ray-tracing technique 
where the effects of scattered emission and re-heating by dust grains are accurately taken into account. To compare 
with observational studies, the dust re-emission spectrum is approximated by a modified black-body function and the 
brightness profile with a Gaussian source. A simplified single-zone model for cores is presented that incorporates on the 
basis of the derived emissivities a first order approximation of their SED. 

Results. Colder dust temperatures are, independent of the core mass, related to a higher pressure both inside and 
around the filament. The pressure-temperature relation for given external pressure is found to be largely independent 
of the true shape of the filament. The calculations show that the radiative transfer leads to a lower emission coefficient 
at 250 /im and to a flatter emissivity law of typically /? < 1.8 in the far-infrared sub-millimeter regime. These effects 
cause the core mass to be underestimated by more than a factor of 2 based on the typical assumptions used in 
observational programs. A larger uncertainty is expected for high pressure regions. Fitting the core using a Gaussian 
source approximation overestimates the flux by ~ 10%. For highly embedded cores and in general for cores in high 
pressure regions, the surface brightness profile is flatter with respect to the profile of the column density and a Gaussian 
profile. These effects can lead to an overestimate of the core size of 10 — 30% based on marginally resolved 250 fira 
observations. 

Key words. Radiative transfer, ISM:dust, extinction, ISM: clouds, Infrared: ISM, Submillimeter: ISM 



1. Introduction 

Condensed cores are compact dense regions of molecu- 
1 lar clouds that have a high column density contrast to 
their surrounding medium. They are regions of star for- 
mation a nd appear to be phy sically very similar to Bok- 
globules (jBok fc Reillvl 1 1 947n . dense, almost spherically 
shaped, highly compact molecular clouds surrounded b y 
• non-molecular thin interstellar gas (jCurrv fc McKcc 2000) . 
The formation of the condensed cores is potentially 
linked to the turbulent motion of the molecular gas that 
produces a self-similar multi-fractal density structure with 
a log-normal distribution of the local density. In this sce- 
nario, the cores are formed out of massive density enhance- 
ments that became self-gravitating and stable against the 
turbulent motion. Turbulence naturally produces a mass 
distribution of clouds similar to that observed in the ISM 
and is understood to be the origin of th e initial mass func- 
tion (I MF) of the stars (see references in lFischera fc Dopital 
(120081) ). 

Because the cores are highly opaque, they are most eas- 
ily studied in the far infrared (FIR) or submm regime by ob- 
serving the dust emission spectrum where the cores are still 
optically thin. Surveys of molecular clouds performed in the 



FIR/submm (as are possible using BLAST, HERSCHEL, 
and ALMA) are essential to our understanding of the very 
early evolutionary process of star formation and will help 
us to unlock the myst e ry of the origin o f the IMF (e.g. 
iNetterfield et all (|2009t ): lOlmi et al.l (|2009H ). The interpre- 
tation, however, of the dust emission spectrum is compli- 
cated by both radiative transfer effects and the unknown 
dust properties in the molecular phase of the ISM. Both 
effects lead to uncertainties in particular in the measure- 
ment of the core masses, which is the main observable de- 
rived by modeling the SED. This causes a major problem in 
gaining information about the efficiency of star formation, 
the stellar mass, or the masses produced within a core of 
a certain mass. The radiative transfer also further affects 
the surface brightness profile of the core as grains in the 
cloud are heated by a gradually changing radiation field 
()Fischera fc DopitalT2008t iNutter et alJl2009ft . 



Observationally, cores display, if described by a single 
temperature model, a large variation in dust temperatures. 
This is a direct consequence not only of their gravitational 
state, but also of their local environment. Some cores are 
exclusively heated by an external radiation field while oth- 
ers have already formed a protostellar object in the center 
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and are therefore also heated from the inside. Cores with 
a warm central ob ject therefore ha ve on average a higher 
dust temperature (|Ohni et al.ll2009|) . 

The SEDs of interstellar clouds and cores hav e been 
extens ively studie d in the past dBernard et al.l 
19931: lEvans et al.l l200fl IStamatellos fc Whitworthl 



1992, 



2003; 



Stamatellos et al.ll2004j : IFischera fc Dopitall2008l) . although 
with different emphases. iBernard et al.l (|1992l) studied the 
radiative transfer using a physical dust model and derived 
the brightness profiles and SEDs of interstellar clouds while 
making simplified assumptions about the scattered light 
and density profiles in clouds (a flat center and a power law 
further out). They concluded that the observed brightness 
profiles cannot be caused by radiative transfer effects alone . 
Evans et al.l (l200lh IStamatellos fc Whitworthl (|2003l) . and 
Stamatellos et al.l (2004) studied the infrared emission from 
cores to obtain more robust estimates of their dust masses. 
Their calculations focused on the radiative transfer as- 
pects and used a simplified dust re-emission model based 
on the mean dust p r operti es that they adopte d from 
lOssenkopf fc Hennins) (| 19941) . IFischera fc Dopital (|2008h 
studied the dependence of the SED derived for stable 
isothermal self-gravitating interstellar clouds on the ra- 
diation field and the external pressure in the interstellar 
medium. These calculations were based on a physical dust 
model of interstellar dust grains. The radiative transfer 
problem was solved by accurately considering the compli- 
cations caused by the sc attered light to provide (as c an als o 
be found in the work of IStamatellos fe Whitworthl ()2003h ) 
profil es of the surface brightn e ss for the whole SED. 

IStamatellos fe Whitworthl (|2003l ) studied the SEDs of 
cores embedded in optical thick dense spherical clouds and 
the present work can be compared to theirs. There are, how- 
ever, basic differences. In their model, the giant molecular 
cloud (GMC) is not modeled as a self-gravitating entity but 
as a homogeneous dense sphere and the pressure inside the 
GMC is assumed to be independent of the extinction of the 
cloud. 

In this paper, I relate the SEDs of passively heated con- 
densed cores to their physical environment. The calcula- 
tions rest upon compar able studies of interstell ar clouds 
that are either spherical (jFischera fe Dopital l2008) or cylin- 
drical (Fischera, in prep.) in shape and are modeled as self- 
gravitating entities. As a first approach, I focus on the ra- 
diative transfer effects and base the SEDs of the cores on 
dust properties in the diffuse ISM. Complications caused by 
the turbulent density structure of molecular clouds are not 
considered and additional studies will be needed to quanti- 
tatively understand how the results presented in this paper 
are affected. The expected effects are discussed. 

To compare the results with observations, the dust re- 
emission spectrum is described by a modified black-body 
spectrum. This provides the corresponding observational 
parameters of the dust emissivity and dust temperature. On 
the basis of these effective values, I present a luminosity- 
temperature relation and a pressure-temperature relation 
of embedded cores. I compare the surface brightness with 
a Gaussian function that is used obser vationally to ex- 
tract the fluxes and the size of the cores (|01mi et al.ll2009l: 
iNetterfield eraHl2009h . 

In Sect. [H I present my model of passively heated em- 
bedded cores, the dust model, and the technique used to 
solve the radiative transfer problem. I also provide for com- 
parison with observational data a single zone model. My 



results of the radiative transfer calculations are presented 
in Sect. 13.11 In Sect. 13.21 I compare my model with the 
modified black-body spectrum and in Sect. 13.31 1 provide 
the surface brightness. In Sect. I discuss the effect of 
my results on uncertainties in the radiation field and dust 
properties, and provide a summary of the main results of 
the paper in Sect. [5] 



2. Model 

2.1. Model of embedded cores 
2.1.1. Model of the cores 

The cores are modeled as isothermal se lf-grav i tating 
spheres , so-c alled Bonnor-Ebert spheres (Ebert 1955; 
lBonnor{IT95fih . which are embedded in extended molecular 
clouds and pressurized by i t s cent ral pressure p c (see refer- 
ences in lFischera fe Dopital (|2008l ^. The cores are naturally 
in different gravitational states. However, for simplicity I 
consider in this work the cores that are in the critical state 
against collapse. 

As described in ICurrv fe McKeel (|2000l ) or 
IFischera fe Dopital (|2008j ). critical stable clouds are 
characterized by a steep density profile at the outskirts 
close to 1/r 2 , that becomes flat at the cloud center. They 
have an overdensity or overpressure p of 14.04 in the cloud 
center where the overpressure p is defined as the ratio 
of the pressure inside the cloud to the external pressure. 
Stable clouds are characterized by a flat density profile and 
lower overpressure, supercritical clouds by a steep density 
profile (also 1/r 2 ) and higher overpressure. Observationally 
the cores are probably in a critical or super-critical state as 
they have a high column density contrast and are therefore 
easy to identify as individual sources. 

The stability of a spherical isothermal self-gravitating 
cloud against gravitational collapse is determined by the 
temperature, the molecular weight, and the pressure out- 
side the cloud. The temperature is considered as an effec- 
tive value that contains, apart from the thermal motion, 
the turbulent motion of the dusty gas. In the critical state 
of the core, its mass is given by 



= 0.688 



T 
10K 



(—) 
V2.36/ 



Pc/k 



2.8 x 10 5 K/cm E 



Mr 



where I have assumed that the pressure p c at the center 
of the GMC is 14 times the mean pressure in the I SM o f 
p/k = 2x 10 4 K/cm 3 proposed bvlCurrv fe McKeel (2000) 
and also adopted bv IFischera fc Dopital ( 20081 ). The cores 
are also assumed to be molecular, which implies a mean 
atomic weight of [i k, 2.36. The cores' radius simply scales 
with the temperature and is given by 



Rc 



0.0345 



T 



10K J V2.36 



Pc/k 



2.8 x 10 5 K/cm 3 



pc.(2) 



In higher pressure regions, cores of the same temperature 
are more compact, have higher density, and are less mas- 
sive. Combining both of these properties provides the mass- 
radius relation for cores 



M r , orfi = 5.78 



Rc 



Pc/k 



0.1 pc/ V 2-8 x 10 5 K/cm 3 



Mr. 



(3) 
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Of prime importance to the radiative transfer calcula- 
tions is the optical depth through the cloud and its density 
profile. The column density through the cloud center for 
a given gravitational state is reflected by the density pro- 
file or overpressure is, independent of the true mass, deter- 
mined by the pressure surrounding the cloud. For a crit- 
ically stable cloud, the column density is, assuming solar 
abundances, given by the expression 



AT H (0) = 5.8 x 10 2 



Pc/k 



2.8 x 10 5 K/cm 3 



(4) 



The column density through the center is difficult ob- 
servationally to measure where the mean column density is 
instead derived. The corresponding mean column density 
of the critical stable core is then given by 



(N H ) = -{N ll ) (0) = 1.6x10 22 



Pc/k 



x 10 5 K/cm 3 



(5) 



where (Nn) (0) is the column density at the center of a 
cloud whose mass is distributed homogeneously inside the 
cloud volume (see Appendix IB. 1[) . 

In higher pressure regions, cores of the same effective 
temperature are less massive, smaller, but more optically 
thick. In higher pressure regions, the dust inside the core is 
heated by a more strongly attenuated radiation field. For 
the same radiation field heating the core, the dust temper- 
atures of the cores in high pressure regions are therefore 
colder than the cores in low pressure regions. 

2.1.2. Model of the filaments 

As the cores, the filaments are modeled as isothermal self- 
gravitating clouds surrounded by a thin interstellar medium 
of pressure p ext . For comparison, both spheres and cylinders 
are considered as approximations of the observed structures 
in the ISM. As self-gravitating entities, they display a den- 
sity and pressure profile. 

It is reasonable to assume that condensed cores form 
out of density enhancements at the centers of the GMC. 
Here, potential cores experience the highest pressures and 
are therefore the locations where clouds can most easily be- 
come self-gravitating. Only those cases are therefore con- 
sidered in the model. For simplicity, the mass and size of 
the core is assumed to be negligible to the mass and the 
extension of the filaments. 

The physical properties of a critical stable spherical 
GMC is described by the same equations as given for the 
cores (Eq. [T] - [5j) . Replacing the pressure by the mean ISM 
pressure of 2 x 10 4 K/cm 3 , the mass radius relation is in 
good a greement with observations of interstellar molecular 
clouds |Larsodfl98l[) . 

To study the whole range of structures from Bok- 
globules to highly embedded cores, I consider for a given 
external pressure different extinction values through the fil- 
aments. Since the central pressure of self-gravitating clouds 
is related to the column density at its center, the central 
pressure is also varied. The relation for spheres and cylin- 
ders embedded in the ISM of our Galaxy is shown in Fig. [T] 
For the same central extinction of a filament, cores situated 
in cylinders are subject to a higher central pressure than 
cores in spheres. 



The range of extinction values considered in the paper 
leads to the situation indicated in Fig. [TJ where the gravi- 
tational state of the spherical filaments varies over a wide 
range from highly stable to, in some cases, highly super- 
critical clouds. For the physical solutions of self-gravitating 
isothermal cylinders, a similar stability classification as for 
spheres does not exist as the cloud can always produce by 
compression a higher gas pressure at the cloud edge. The 
cylindrical filaments can therefore for all cases studied in 
the paper be regarded as stable. 

Cylinders have a maximum mass-line density given by 
M cy \/L = 2K/G, where M cy \ is the cloud mass, L is 
the cloud length, and G is the gravitational constant. 
The constant K, given by k^T/{p,m u ), is the one gener- 
ally introduced in isothermal cloud models where /cb is 
the Boltzmann constant, m n the atomic mass unit, and fi 
the mean molecular weight. This mass-line density is only 
reached in the limit of infinite overpressure (Fischera, in 
prep.). By construction, the mass-line density of cylinders 
with finite overpressure is lower than the asymptotic value. 
If we consider two clouds with the same mass-line density, 
a higher overpressure corresponds to a lower (effective) gas 
temperature. A change in the external pressure clearly has 
no consequence on the overpressure. 

The situation is fundamentally different for spherical 
clouds where the maximum mass is identical to the critical 
stable mass given by M sp h oc {T/ p) 2 / ,/Pext, where a higher 
external pressure will change the gravitational state and 
can make the cloud unstable. While a cooling of a spherical 
cloud will consequently lead to a critical and finally unsta- 
ble cloud it will lead only to a higher gravitational but still 
stable state in the case of cylinders. I note that - although a 
supercritical spherical cloud can be constructed - this state 
can not be obtained by means of simple cooling. If pressure 
equilibrium between cloud and the surrounding medium 
were assumed cooling of a critical stable cloud would lead 
to a decrease in the pressure at the outskirts. A spherical 
super-critical cloud assumed to be in pressure equilibrium 
with the surrounding medium also corresponds to a cloud 
temperature above the critical temperature (Fischera, in 
prep.). 

Cylindrical self-gravitating clouds of high overpressure 
are characteri zed by a steep density profile at the outskirts 
with p oc r~ 4 ()Stod61kiewiczlll963t[Ostrikerlll964D . The ob- 
ser yed density profi l e of int erstellar clouds as pointed out 
by iFiege fc Pudritd (|2000al) appears to be less steep and 
more closely represented by a profile p oc r~ 2 . Instabilities 
m ay explain these r e sults as e.g. discussed bv lLarsonl (£l985) 
or lFiege fc Pudrita (|200 0b). which prevent the clouds pro- 
ducing high overpressures. However, cases with high over- 
pressure are still useful in the framework of the radiative 
transfer problem and are therefore included in the model. 

The column density at the center of self-gravitating 
cloud models is given by 



JVh(O) = 



Pc 



4ttG 



(6) 



where w is the potential, G the gravitational constant, 
Pi the weight of element i in units of atomic mass units 
m u , and z c \ is determined by the pressure equilibrium be- 
tween the cloud and the surrounding medium p{z c \) = 



p c e~ 



p ex t. The column density at the center of a 
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Fig. 1. Column density at the center of a spherical (black) 
and cylindrical (grey) isothermal self-gravitating cloud as 
a function of the overpressure p c /Poxt, where p c is the cen- 
tral pressure and p C xt is the external pressure. The dashed 
lines indicates the asymptotic behavior for high overpres- 
sure where Nn(0) oc ^fp~ c . The external pressure is assumed 
to be consistent with the mean pr essure of 2 x 10 4 K/cm 3 
in the ISM of the Milky Way (|Boulares fc Coxl [1991 . 
The grey dotted lines show the values for a critical sta- 
ble sphere with an overpressure of 14.045. The diamonds 
mark the values considered in the paper and correspond 
to Ay = [1,2,4,8,16], where the extinction Ay is de- 
rived from the column density using the d ust-to-gas ra- 
tio E (B - V)/N K = 5.8 x lCT 21 mag/cm- 2 (jBohlin et al.l 
Il978f) and the ab solute to relative extinction Ry = 3.1 
(|Fitzpatricklll999h . 

core given in Eq. 0] is a special case. In the limit of high 
overpressure, the integral in Eq. [6] becomes a constant and 
the column density solely a function of the central pressure 
p c . The asymptotic behavior is given by 



iV H (0) ~ 1.5 x 10 21 £. 
where 



Pc/k 



2 x 10 4 K/cm E 



cm 



3.0279 
2.2214 



(spheres), 
(cylinders) . 



(7) 



(8) 



Because of this characteristic, in the case of high overpres- 
sure, the radiative transfer problem and as a consequence 
the dust temperature of the cores become largely indepen- 
dent of the external pressure. 

2.2. Dust model 

The dust is assumed to be a specific composition of 
graphites, silicates, iron grains, and PAH molecules where 
each dust component has a give n size distribution. I use th e 
same dust model parameters as lFischera fc Dopital ((2008), 
which were chosen to reproduce the observed dust proper- 
ties of the diff use phase of the IS M including the mean ex- 
tinction curve (|Fitzpatricklll999D . the diffuse dust emissio n 
derived with the DIRBE experiment (jArendt et al.lll998h . 
and the depletion of key elements from the gas phase. 

The dust re-emission spectrum and the PAH emission 
spectrum are derived by taking into account the stochastic 



nature of the heating process for individual grains. This 
heating process is of particular importance for small dust 
particles where the thermal energy is typically lower than 
the absorbed photon energy, leading to large temperature 
variations. 

There are some di fferences between my model and that 
of Draine fc Lil (|2007h as mentioned in iFischera fc Dopital 
(12008ft . First, the dust albedo in the optical is lower; sec- 
ondly, the emission at shorter wavelengths where the dust 
emission spectrum is supposedly dominated by stochastic 
emission is enhanced by additional iron grains; and thirdly, 
the graphite grains are restricted to relatively large sizes 
with radii larger than 0.01 /im to avoid any contribution 
to the extinction bump at 2200 A, which is entirely ex- 
plained by PAH molecules. The dust-to-gas ratio in mass 
is r\ — 0.59%. As the grains possibly accrete gas atoms or 
molecules, the dust-to-gas ratio might indeed be slightly 
higher. 



2.2.1. Mean extinction coefficient at long wavelengths 

At far infrared wavelengths, where most of the grain cooling 
occurs, the extinction is close to the absorption probability. 
The mean cross-section for extinction per hydrogen atom 
in the model is given by the sum over the different dust 
components and the integrations over the individual size 
distributions 



(cr) /h 



da/, (a)™ 2 Q c ^(i,a), 



(9) 



where <5^ xt (i,a) is the extinction coefficient at wavelength 
A of a dust grain of radius a and composition i. The size 
distributions ,fi{a) are normalized such that Ci — ni/n>H is 
the relative amount of a key element th a t is pa rt of the dust 
composition i (see IFischera & Dopital (|2008D 1. The same 
expression as Eq.|9]is used to determine the scattering cross- 
section and the absorption cross-section. The dust albedo 
is given by the ratio uj\ — (C| ca ) / (CJ**). 

The mean cross-section for the dust extinction assumed 
in the model calculations is shown in Fig. [5] The values in 
the FIR/submm re gime are close, albeit slightly lo wer, than 
the ones derived bv IWeineartner fc Drainel pOOl). The be- 
havior is predominantly determined by graphite, PAH, and 
silicate grains. The relative importance of the components 
changes as a function of wavelength. While the extinction 
in the optical is caused mainly by carbonaceous grains, sil- 
icate is the dominant absorber from 10 /im up to 60 /im. In 
the FIR/submm regime, the extinction is caused almost in 
equal parts by silicate and graphite grains. The extinction 
curve at wavelengths above 40 /im, in the spectral region 
where the bulk of the cold dust emission occurs, is approx- 
imated well by a power law. The extinction coefficient per 
unit mass is given by 



(Pd) /H ' 

where 

(p d ) /H « 1.384 x 10~ 26 



H — atom 



(10) 



(11) 



is the mean dust density per hydrogen atom (Fischera, 
in prep.). The long- wavelength behavior above 60 /im to 
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Fig. 2. Mean cross-section of the dust extinction per hy- 
drogen atom as assumed in the radiative transfer model. Its 
functional dependence on wavelength is compa red with the 
obser ved mean extinction for the diffuse ISM ([Fitzpatrickl 
Il999f) (although extend ed to the FIR) and the mean extinc- 
tion curve provided by IWeingartner fc Draind ([2001[ ) (cor- 
rected by a factor 0.93). Also shown are the contributions 
of the several dust components in the model. 

1000 /im is (within 4% accuracy to the dust model) de- 
scribed by 

/ \ \- 2 024 rm 2 
« 5.232 ( — ) — . (12) 

The power-law behavior of the dust model also extends to 
longer wavelengths. The approximation up to 10 4 fxm is 
accurate to within 10%. 

The extinction coefficient per unit length is fc| xt = 
({Cf*)/H)nu. 

2.3. Radiative transfer 

The filaments are assumed to be heated by a n isotropic in- 
terste llar radiation field (ISRF) adopted from Mathis ct al. 
(| 1982T) . The spectrum is characterized by a sharp cutoff at 
13.6 eV and contains, apart from the stellar component, no 
re-emission from dust grains. 

The radiative transfer is solved by assuming that the 
mass and size of the core is negligible compared to the 
mass and size of the filaments, the cores are embedded in. 
Under these circumstances the radiative transfer problem 
can be solved in two separate steps: The radiative transfer 
through the filament and the radiative transfer through the 
core. The centrally located cores are heated by the mean 
intensity J\ (0) at the center of the filaments derived for the 
filaments without cores. For simplicity, the radiation field 
caused by the attenuated ISRF, the scattered emission, and 
the thermal dust emission from the filaments is assumed to 
be isotropic. 

For fixed dust properties inside the filaments and cores, 
the radiative transfer problem (and therefore the SEDs of 
the cores) is determined by the pressure of the gas sur- 
rounding the filaments and the overpressure p c /p C xt in the 
filaments. The shape of the filaments has an additional ef- 
fect on the radiative transfer because it determines the ra- 
dial density profile, the extinction- central pressure relation 



fSect. l2~l~2")) and the attenuation of the radiation inside the 
filaments. The influence of the shape of the filament on 
the SEDs of the cores is explored by considering both self- 
gravitating spheres and cylinders. 

2.3.1. Accurate solution of the SEDs of the cores 

For the given simplifications, the radiative transfer prob- 
lem as a function of external pressure and overpressure 
is accurately solved usin g a ray-tracing technique (see 
iFischera fc Dopital ([200 8)). This technique provides the 
angle-dependent radiation field inside the cloud caused by 
the attenuated radiation, the scattered radiation, and the 
thermal emission from dust grains. In the procedure, the 
cloud is divided into a number of optically thin shells. The 
impact parameters of the individual rays, which is their 
nearest distance to the cloud center, is given by the mean 
distance to the shells. For cylindrical geometric calcula- 
tions, rays at 30 additional angles around the main rotation 
axis are considered (Fischera, in prep.). Effects caused by 
non-isotropic scattering and multiple scattering events are 
accurately included in the calculation as well as the effects 
of the re-heating of the dust grains and subsequent emis- 
sion. The polarization of the scattered emission is assumed 
to have a negligible effect on determining the SED and is 
therefore ignored. 

2.3.2. Approximate solution for the SEDs of embedded cores 

When observations are analyzed, the dust re-emission spec- 
trum of condensed cores is typically modeled by assuming 
a single dust temperature where the dust re-emission spec- 
trum is approximated by a modified black-body spectrum. 
The SED (in units W/m 2 /A) from the core of dust mass 
Afdust at distance D corc is given by 

5 , ust= M , tKx Bx{Tdnsth (13) 
core 

where Td us t is the dust temperature, -BA(2dust) the Planck- 
function, and the emission coefficient per unit dust 
mass. The emissivity is commonly approximated by a power 
law k™ 1 = Ko(A/Ao) _/3 by making certain assumptions 
about the dust emissivity at wavelength An an d the 
power /? dNetterfield et all [20091 lOlmi et al.l [2009). The 
model is equivalent theoretically to a single-zone model for 
the dust temperature in the core. However, as shown in this 
paper, to obtain the correct shape of the dust re-emission 
spectrum and the correct dust mass, the intrinsic values 
of the dust emission coefficient need to be replaced by the 
appropriate effective values to take into account not only 
the variation in the emission spectrum with both size and 
composition but also the radial variation in the grain tem- 
peratures. I use this approximation to obtain for the given 
dust mass and theoretical SED the appropriate effective 
parameters. 

The single-zone model allows us not only to interpret 
the observations but also to test efficiently a wider range 
of model parameters. It is therefore a more practical model 
to use than time-consuming ray-tracing calculations. The 
model is described in detail in Sect. [2J In the approximate 
solution of the radiative transfer problem based on energy 
conservation the effects of both scattering and re- heating by 
dust grains inside the filaments and the cores are taken into 
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account. However, in the considered cases the contribution 
of self-absorption to the total heating rate inside the clouds 
is negligible. Of greater importance is the contribution of 
dust emission from the filaments to the total heating of 
highly embedded cores. In the model approximation this 
additional heating is provided as a lower limit. 

3. Results 

The dependence of the effective parameters is analyzed for 
three assumptions about the external pressure. Apart from 
the mean pressure in the ISM (2 x 10 4 K/cm 3 ), situations 
where the filaments are embedded in a higher pressure gas 
(10 5 and 10 6 K/cm 3 ) are also considered. For each consid- 
ered value of external pressure, the SEDs of the cores are 
derived for a range of extinction values Ay at the center 
of the filaments. For both spherical and cylindrical geome- 
tries for the filaments, I considered the extinction values 0, 
I, 2, 4, 8, and 16. For cores in spherical filaments, I also 
assumed higher extinction values (32 and 64 mag). In these 
calculations, the energy is conserved in all cases to better 
than 0.2%. 

In the following section 13.11 I describe the individual 
steps of the radiative transfer calculations. The results ob- 
tained using a modified Planck-function are described in 
Sect. [321 

3.1. The theoretical dust emission spectrum 

3.1.1. The mean intensity at the filaments' center 

The first result obtained using the ray-tracing technique is 
the mean intensity at the center of the filaments. The rela- 
tive contribution of the attenuated external radiation field, 
the scattered emission, and the re-emission from dust grains 
located in the filaments to the SED is shown in Fig. [3] for 
two extinction and pressure values. For the ISM pressure 
(Pext/k — 2 x 10 4 K/cm 3 ), the extinction value Ay = 8 mag 
in the spherical filament is close to the value of a criti- 
cal stable sphere with overpressure 13.2, while the higher 
extinction value, Ay = 64 mag, corresponds to a super- 
critical spherical cloud with high overpressure (68.25). In 
both cases, the clouds are characterized by a steep den- 
sity profile in their outskirts. In the high pressure region 
(Pext/k = 10 6 K/cm 3 ), the high extinction value is close to 
the corresponding critical value (Ay ~ 59 mag), while the 
low extinction value corresponds to a sub-critical spherical 
cloud. The clouds in the high pressure region are charac- 
terized by flat or steep density profiles with low and high 
overpressures (1.30 and 16.38, respectively) at the cloud 
center (see tabic [1]) . 

The SED of the center of the filaments is primarily de- 
termined by the radial extinction and in parts by the geom- 
etry. The effects caused by the radial density profile play 
only a minor role in particular for the dust re-emission spec- 
trum. The properties are of importance for interpreting the 
dust re-emission from cores and are described in the follow- 
ing. 

Because of simple geometric effects, the radiation is, for 
identical values of Ay, more strongly attenuated by cylin- 
drical than spherical geometry. The strength of the scat- 
tered light that dominates the SED below its maximum 
depends to a certain degree on the shape of the filaments. 
For the considered high extinction values, most of the scat- 



tered emission in the UV and optical is generated in a skin 
layer of the filament, so that the geometric effects produce 
a lower amount of scattered emission at the center of the 
cylinders. The situation is different at longer wavelengths 
where the clouds become optically thin and the probabil- 
ity of scattering events is higher in cylinders than spheres. 
Despite its low intensity for the cases shown in Fig. [31 the 
scattered light is responsible for appro ximately half of the 
dust heating in the filaments' center ([Fischera Sz Dopital 
I2008h . 

At the high extinction values considered, the amount 
of PAH emission produced in the filaments is rather in- 
sensitive to the extinction value since most of the flux re- 
sponsible for heating the molecules is absorbed. The emis- 
sion originates mainly in a layer close to the surface where 
the UV and optical light most responsible for heating the 
molecules is strong. The radiation in this layer is more in- 
tense in low pressure regions causing for a given extinction 
a slightly stronger PAH emission at the filaments' center. 
In the case of Ay — 64 mag, the PAH emission itself ex- 
periences a certain extinction that is most evident in the 
larger minimum at 10 /im, which is caused by silicate ab- 
sorption. It can be seen overall, that the ratio of the PAH 
emission relative to the dust emission peak is considerably 
lower than in the diffuse ISM (Fig. [15]). 

The shape of the dust re-emission in the FIR/submm- 
regime depends for a given extinction and external pressure 
on the shape of the filaments. Radiation from cold dust 
grains located in the interior of the filament is higher in the 
case of a cylinder because of the infinite column density 
along the major axis. At relatively low extinction values 
(Ay ~ 1 mag) where the UV and optical light can pene- 
trate deeper into the cloud, this also affects the intensity of 
the PAH emission and the warm emission of stochastically 
heated small grains. The dust emission peak at higher ex- 
tinction values also appears to be broader because of the 
larger temperature variations and colder dust temperatures 
at the filaments' center. 

3.1.2. Grain heating 

The dust heating 

h(r) = J d\(kf s /n n )^J x (r) (14) 

at the core caused by the mean intensity in the filament 
varies strongly inside the core (as seen in Fig. @| because 
the light, especially in the UV and optical, is strongly atten- 
uated by dust scattering and dust absorption towards the 
core center. This radial variation is even stronger in the 
high pressure region (p e *t/k = 10 6 K/cm 3 ), where cores 
experience higher pressure and are more opaque. 

In Fig. 21 the heating is divided into several components 
corresponding to the various radiation fields inside the core. 
For cores in filaments with Ay = 8 mag, the heating is dom- 
inated almost in equal parts by the attenuated ISRF and 
the scattered light. As visible in the figure, inside the core 
the scattered radiation field is mostly caused by photons 
scattered inside the core. Deep inside highly opaque cores 
and in the case of highly embedded cores the scattered con- 
tribution decreases relative to the attenuated ISRF as the 
heating shifts to longer wavelengths where the dust albedo 
is lower. 
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Fig. 3. Mean intensity at the center of spherical (black solid lines) and cylindrical (grey solid lines) filaments heated by 
the ISRF (long dashed lines). The filaments are embedded in a medium withp cxt /fc = 2 x 10 4 K/cm 3 (top) and 10 6 K/cm 3 
(bottom). The extinction through the filaments' center is either Ay = 8 mag (left-hand figure) or Ay — 64 mag (right- 
hand figure). Shown arc the different components of the mean intensity: scattered emission (dotted lines), attenuated 
light (short dashed lines) , and dust re-emission (dashed-dotted lines) . The emission peak in the millimeter regime is the 
CMB radiation. 



In case of a strongly attenuated external radiation field, 
the heating by dust emission from the filament and the 
core becomes an important component. This applies ei- 
ther for the central region of highly opaque cores and for 
highly embedded cores in general. The size of the contri- 
bution in reality depends certainly on the dust properties 
in the dense phase of the ISM. The heating derived here 
may be an upper limit since most of the heating is caused 
by the PAH emission which is possibly absent in realistic 
clouds. The heating by self- absorption plays a larger role 
in t he case of non-embedded cores in high pressure regions 
(see lFischera fc Dopital (|2008D V For the considered cases in 
Fig. |4] of embedded cores, the contribution is smaller than 
~ 10%. Far more important is the heating caused by the 
dust emission from the filaments. At Ay = 8 mag, this 
component produces at the center of the cores embedded 
in spherical filaments between 13% and 25% of the total 
heating rate. The contribution is even higher if the cores 
are instead embedded in cylinders (20 — 35%). For highly 
embedded cores, the dust in the whole interior of the cores 
is predominantly heated by the dust emission from the fil- 
aments. This additional heating limits the variation in the 



heating rate inside the core. I note that the heating com- 
ponent is higher in the low pressure region. 

The effect of the external pressure and the extinction 
of the filament on the temperature variation in the heat- 
ing rate inside the core is further visualized in Fig. [5] where 
the heating rate in the center is compared with the heating 
rate at the cores' edge. For comparison with simplified ra- 
diative transfer models, the heating rate is also related to 
a dust temperature. The heating rates at the cores' edge 
are determined in a preliminary way for a given filament 
shape by the extinction Ay through the filaments' center. 
Since the cores are more compact for a given extinction 
in higher pressure regions that produce a lower intensity 
at the cores' edge, the heating rate decreases slightly to 
higher pressures. The effect reaches a maximum of only a 
factor of two. Clearly visible in the figure is the dependence 
of the overpressure on the extinction of the filaments for 
a given external pressure. While the overpressure changes 
drastically in the low pressure regions also in case of low 
extinction, the effect is small for filaments in high pressure 
regions. 
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Fig. 4. Dust heating inside cores embedded in spherical (black curves) or cylindrical (grey curves) filaments. The visible 
extinction Ay at the center of the filaments is either 8 mag (left-hand figure) or 64 mag (right-hand figure). The filaments 
are pressurized either by the mean pressure in the ISM (top) or by a medium with p ^t/k = 10 6 K/cm 3 (bottom). The 
total heating rates (solid curves) are separated into different components corresponding to the different radiation fields 
in the core: Attenuated ISRF (long dashed curves), total scattered light (dotted curves), light scattered in the core 
(short dashed curves), dust emission from the filament (dashed-three dotted curves), and dust emission from the core 
(dashed-single dotted curves). The heating by dust emission produced inside the core is multiplied by a factor 5. 



The lowest dust temperature in the core is strongly de- 
pendent on the pressure surrounding the core, which de- 
termines its extinction. However, there are some noticeable 
systematic effects caused by the shape of the filament and 
the external pressure. At a certain pressure p c , the mini- 
mum temperatures are naturally lower for embedded cores 
in low pressure regions p cxt than for non-embedded cores in 
high pressure regions. The behavior at high pressure values 
p c for a given external pressure reflects the differences in 
the dust emission inside the filament noticeable in Fig. [4] 
For given external and internal pressures, the dust at the 
center of cores is for example slightly warmer, if they are 
embedded in cylindrical filaments as the intensity of the 
dust emission is higher in cylinders than in spheres. 

As seen in Fig.HJ the qualitative behavior of the relative 
variation of the heating rate or the temperature inside the 
core depends on the pressure surrounding the filaments. In 
high pressure regions, the variation in the dust temperature 
inside the cores decreases if the core becomes more embed- 
ded and shielded from the ISRF. As the cores become more 
optically thin to the reddened radiation, the radial varia- 



tion in the heating rate and consequently that of the dust 
temperature become smaller. 

The situation in the low pressure region is more com- 
plex as the increase in the extinction of the filament is ac- 
companied by a higher pressure at its center. The embed- 
ded cores are not only heated by a reddened radiation field 
but are also more opaque than non-embedded cores. Above 
Ay = 1 mag, this leads, assuming a mean ISM pressure, 
to an increase in the dust temperature variation with cen- 
tral pressure. As will be shown this causes a flattening of 
the overall SED with extinction (Sect. 13.2. ID . The varia- 
tion in the dust temperature inside embedded cores is not 
necessarily lower than for non-embedded cores. Based on 
my calculation, the cores in molecular clouds of typical ex- 
tinction Ay = 8 display even slightly larger temperature 
variations than non-embedded cores. 

Independent of the assumed external pressure, the vari- 
ation in the heating rate decreases for highly embedded 
cores within filaments with Ay 3> 16 mag, where the inte- 
rior of the cores becomes predominantly heated by the dust 
emission from the filaments. 
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3.1.3. Grain temperatures 

In the model, the SEDs of the cores are produced by the 
re-emission of individual grains of different sizes and com- 
positions that have unique temperature and emission be- 



haviors. The temperature of the individual grains as a func- 
tion of radius is shown in Fig. \E[ The grain sizes are large 
enough for their temperature distribution to be close to a 
narrow distribution around the equilibrium temperature. 
The radial variation in the temperature is most visible for 
the smallest grains, which are heated predominantly by the 
strongly diminished UV and optical light. The figure also 
shows the dependence of the dust temperature on compo- 
sition. Small iron grains are considerably warmer than sil- 
icate or graphite grains because of the high absorption in 
the UV a nd optical and the low emission probability in the 
infrared (Fis cheral 12004 ). In cores embedded in filaments 
with Ay — 8 mag that are pressurized by the mean ISM 
pressure, the temperatures of large grains lie between 8 K 
and 15 K. I note that the temperature variation for dif- 
ferent grain types is smaller in both the core interior and 
more opaque filaments where the grains are heated by a 
more reddened radiation field. The physical explanation is 
the similar absorption and emission behavior of the grains 
at long wavelengths, which for most grains is closely de- 
scribed by a power law with Q^ bs / cm ^ a /X 2 where a is 
the grain size. In this limit, the ratio of the cooling to the 
heating rates and therefore the grain temperatures become 
the same. 

3.1.4. The SEDs from cores 

Figure [7| shows the obtained SEDs of cores, which are 
embedded in spherical filaments with Ay = 8 mag and 
Ay = 64 mag. The SEDs of cores embedded in cylindri- 
cal filaments have the same properties but slightly lower 
PAH emission and lower dust temperatures for the same 
extinction values. 

The SEDs reflect the dependence of the opaqueness of 
the cores on both extinction and external pressure as dis- 
cussed in the previous sections. For highly embedded cores, 
the SEDs depend to most parts on the extinction of the fil- 
ament and is largely independent of the external pressure 
(Sect. I2.1.2p . The cores are highly opaque and are heated 
almost in equal part by the attenuated stellar radiation field 
and the PAH emission from the filaments (see Appendix \K[ 
Fig. IA.l]l . The minimum around 10 /im is enhanced by sil- 
icate absorption. 

At the blue side of the spectrum, the emission is dom- 
inated by light scattered within the core. For a given ex- 
tinction Ay, the scattered radiation is only mildly depen- 
dent on the pressure since most of the scattered light origi- 
nates in a layer close to the surface of the core (see Fig. 0]). 
Multiple scattered light only becomes important at wave- 
lengths A < 3 /im. 

The bulk of the dust re-emission is closely described by a 
modified black-body function. At shorter wavelengths, the 
spectral shape is dominated by the emission from stochas- 
tic dust emission, small grains, and PAH molecules. The 
effect of the dust attenuation inside the core is clearly vis- 
ible in the PAH emission. Because of the low intensity of 
the UV/optical, the ratio of the PAH emission to the dust 
emission peak is lower than in the diffuse ISM or inside the 
filaments (Sect. [3~TTT|) and decreases systematically with ex- 
tinction (Fig. [5]). For highly embedded cores, the reddened 
radiation affects the spectral shape of the PAH emission 
spectrum as the photons are not sufficiently energetic to 
heat the molecules to high temperatures. 
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Fig. 6. Dust grain temperatures for a number of different sizes and dust compositions of cores embedded in spherical 
(black lines) or cylindrical (grey lines) filaments. The iron temperature is divided by a factor of 2. The effective values 
for a single temperature model of the dust grains is shown for comparison (solid curves) where the dust emissivity is 
approximated by a simple power law fEq. I12I) . The external pressure p eyi t/k outside the filaments is assumed to be either 
2 x f0 4 (top) or 10 6 K/cm 3 (bottom). The extinction through the filaments' center is either 8 mag (left-hand figure) or 
64 mag (right-hand figure). 



For a given extinction Ay , the PAH emission of embed- 
ded cores is not strongly dependent on the external pressure 
as the heating flux for the PAH molecules is largely ab- 
sorbed in the core. The dust emission, on the other hand, 
appears to be colder and stronger in higher pressure re- 
gions as the grains in the core's center are heated by a 
more strongly attenuated radiation field with a larger frac- 
tion of the radiation being absorbed inside the core. There 
is, however, a limit to the total luminosity given by the 
ISRF heating the filaments (see Sect. I3.2.2I) 

The effect of the extinction for a given external pres- 
sure Pext/k on the dust emission from cores is shown in 
Fig. [51 which covers the range from non-embedded (Bok- 
globules) to highly embedded cores. In low pressure regions, 
a higher extinction of the filaments leads to an enhanced 
emission in the Rayleigh- Jeans limit of the spectrum. This 
behavior is responsible for the pressure-extinction relation 
of the filaments. As I have shown (Sect. I57L2"]) . this relation 
significantly decreases the heating rate inside the core and 
produces above Ay ~ 1 mag a larger temperature variation 
inside the core. The effect disappears if the extinction does 
not lead to a strong increase in the central pressure as in 



the high pressure region below Ay ~ 16 mag. The lower 
amount of radiation from embedded cores is then largely 
compensated by the shift to colder dust temperatures. 

3.2. The effective parameters for the SED 

For direct comparisons with observations, the theoretical 
dust emission spectrum at long wavelengths is fitted by a 
modified black-body spectrum having Tdust, ft, and kq as 
free parameters. The fit is achieved using a non-linear x 2- fit 
with 

2/ q rr \ (F\t ~ Sx^Ko, (3,T dust )) 2 , . 

X {Ko,P,T dust ) = > -2 , (15) 

where F\ t and S\ i refer to the mean theoretical flux and 
the mean fitted flux of the modified black body spectrum 
(Eq. [T2]) over wavelength bins i, respectively. 

As can be seen in Fig. [71 the dust emission spectrum at 
shorter wavelengths from the maximum deviates consider- 
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Fig. 7. SEDs of cores located at the center of a spherical filament with a total central extinction Ay = 8 mag (left-hand 
figure) and Ay = 64 mag (right-hand figure). The filaments are embedded in a medium with p ex t/k = 2 x 10 4 K/cm 3 
(top) or 10 6 K/cm 3 (bottome) and are heated by an unattenuated interstellar radiation field (grey solid line). The derived 
dust re-emission spectrum is shown as a solid black line, the radiation heating the core as a long dashed line, and the 
attenuated emission leaving the core as a dashed-dotted line. Also given is the emission produced by photons scattered 
inside the core (black dotted lines). The grey dotted lines show the individual components of multiply scattered photons. 
The top curve shows singly scattered emission, the curves below emission scattered twice and so on. The dust emission 
spectrum is fitted to a modified black-body function (short dashed line). 
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Fig. 8. SEDs of condensed c ores embedded in the center of spherical self-gravitating filaments illuminated by the ISRF 
given bvlMathis~et al. (1983). The external pressure p ex t/k is either 2 x 10 4 (left-hand figure) or 10 6 K/cm 3 (right-hand 
figure). The SEDs are labelled with the assumed extinction Ay through the filaments' center. 
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ably from the simple modified black-body spectrum. The 
fit is therefore restricted to the interval 

[1400/T dust [#],10 000] fim. (16) 

The derived parameters for the modified black-body ap- 
proximation are listed in Table [TJ where I also provide 
the luminosity of the modified black-body function and 
the integrated luminosity. In the model, at low extinction 
values, a larger fraction of the light is absorbed by PAH 
molecules and small dust particles and subsequently re- 
emitted at shorter wavelengths. This causes the modified 
black-body function to underestimate the total dust emis- 
sion. The PAH/small dust particle contribution disappears 
gradually towards highly embedded cores as the UV and 
optical light becomes almost entirely absorbed inside the 
filaments heating the small dust particles. The last two 
columns give the FWHM and the total integrated flux of 
an Gaussian source approximation for the brightness profile 
(see Sect. OS]). For ft, k , L/M COIe , and the FWHM, I pro- 
vide polynomial fits as functions of the dust temperature 
(Table [5]). I emphasize that the derived relations depend 
on the specific dust properties, and to an extent, on the 
radiation field. The relations provided are appropriate for 
dust with the properties of diffuse ISM dust. 

As an estimate of the accuracy of the achieved fit, I list 

i/l/Xred where Xrcd = xV^W and A frcc = N ~ 3 is the 
number of free parameters. The accuracy is within 2.5% to 
4% with the largest deviation occuring around the peak of 
the dust re-emission spectrum. The fit becomes worse for 
higher extinction values because of the larger temperature 
variations in the core. Closer agreement is achieved again 
at extinction values larger than Ay = 64 mag. Because 
of the larger temperature variations inside the core for a 
given central extinction of the filament, the agreement is 
also less good in high pressure regions. I note that in part 
the disagreement is also caused by the optical properties 
used in the radiative transfer calculations, whose properties 
differ, as those of graphite grains, strongly from a simple 
power-law behavior in the FIR. 

In low pressure regions, the filaments and the condensed 
cores become less optically thin. In the limit of a very low 
pressure region, the dust emission from the cores is identical 
to the one in the diffuse ISM. To obtain the corresponding 
effective values of the modified black-body approximation 
for the diffuse dust emission, I considered a spherical self- 
gravitating cloud with Ay = 0.01 mag (assuming p/k = 
2 x 10 4 K/cm 3 ). 

The fit of the dust re-emission spectrum is not as good 
as for emission from the dense cores because the contri- 
bution of stochastically heated small grains at the blue 
side of the emission spectrum is stronger. The derived ef- 
fective temperature is also at Td us t = 19.86 K relatively 
high. The parameters of the emissivity are ft = 1.83 and 
n = 3.40 cm 2 /g. 

If the grains are heated by the non-attenuated ISRF, 
the ratio of the luminosity of the dust emission to the total 
mass of gas and dust is close to '1' (in solar units). The ratio 
obtained with the modified black-body fit is 0.736 Lq/Mq. 
For the total dust emission spectrum, the ratio is higher at 

itotal/Mcoro = 1-149 Lq/Mq. 

The properties of the effective emissivity, the luminosity, 
and the temperature are described in the following subsec- 
tions. It is found overall that the dependence of the effective 
parameters on dust temperature is not strongly affected by 



the shape of the filaments. For the same overpressure, the 
central extinction of a cylinder is indeed lower. However, 
when averaging over the orientations the different geome- 
try produces in the mean an extinction close to the one 
of a spherical filament. The dust temperature for a given 
overpressure is therefore closely the same. 

As shown in Sect. 12.1.21 at high overpressure, as for 
highly embedded cores, the radiative transfer problem does 
not strongly depend on the external pressure. Therefore, 
the effect of the external pressure on the parameters of the 
modified black-body spectrum weakens towards colder dust 
temperatures from cores. 

3.2.1. The effective dust emissivity 

In Fig. 121 I present the spectral index ft and the emissivity 
kq at Ao = 250 /Ltm as functions of the dust temperature 
of the cores. In all cases, the emissivity is lower than the 
corresponding mean value. Furthermore, the emissivity is 
less steep with ft < 1.8 in all cases. The parameters ft and 
kq are lower in higher pressure regions. The origin of the 
flatter emissivity is discussed in Sect. 14.11 

As the figure shows, the emissivity varies as a function 
of the dust temperature. This dependence decreases with 
the external pressure. In the high pressure region, the emis- 
sivity for Ay < 16 mag shows only a mild dependence on 
dust temperature with ft w 1.5. In lower pressure regions, 
the emissivity is flatter in colder cores. For the ISM pres- 
sure, ft decreases from 1.75 to 1.55. For typical GMCs with 
Ay = 8 mag, ft is close to 1.6. The flattening reflects the in- 
creasing optical depth with extinction Ay through the fila- 
ments as shown in Sect. 12.1121 This dependence is strongest 
in the lowest pressure region. For highly embedded cores 
with Ay > 32, the apparent emissivity steepens again and 
seems to approach that corresponding to the intrinsic dust 
properties. The behavior is discussed in Sect. 14.11 

As can be seen from their functional dependence on dust 
temperature, the values of ft and kq are strongly correlated. 
A low ft is accompanied by a low Kq as shown in Fig. 1101 
This correlation is most apparent for regions with p cx t/k < 
10 6 K/cm 3 and extinction values Ay < 16 mag, where 
the data points exhibit only very small deviations from the 
fitted polynomial functions. 

3.2.2. Luminosity of passively heated cores and the 
P - T-relation 

The model relates the L/M-ratio and the effective dust 
temperature of the core to the central pressure in the fil- 
aments surrounding the core (Fig. QT]). For a given shape, 
ISRF, external pressure, and extinction, the ratio is inde- 
pendent of the true size of the core since both the luminos- 
ity and the mass, by mean of the mass-radius relation, are 
proportional to R^ mc - 

If the flux heating a non-embedded core is completely 
absorbed, the luminosity becomes L — 47ri?^ orc it J ISRF , 
where J ISRF is the integrated flux of the ISRF. Ignoring 
the emission of t he CMB, the ISRF provided by 
iMathis et al.l (|1982f) has an energy density u = ^-J — 
0.46 eV/cm 3 . Adopting the mass-radius relation (Eq. [3]), 
the luminosity-mass ratio is given by 
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Fig. 9. j3 and kq as functions of the effective dust temperature T^ us t for the dust re-emission spectrum for cores assuming 
diffuse dust. The cores are either embedded in spherical or cylindrical filaments. The external pressure p ex t/k outside the 
filaments is assumed to be 2 x 10 4 , 10 5 , or 10 5 K/cm 3 . The corresponding theoretical values labelled by the considered 
central extinction Ay of the filaments are shown either as diamonds, triangles, or squares. The values obtained for a 
given external pressure are fitted with a polynomial function shown as solid black lines. 
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Fig. 11. L/Af -ratios (left-hand figure) and the effective dust temperatures (right-hand figure) for embedded passively 
heated cores. The values obtained using ray-tracing (symbols) are compared with the approximate solutions for cores 
in spheres (solid curves) and cylinders (dashed curves). The different symbols refer to the external pressures (shown as 
vertical grey lines) assumed in the calculations: diamonds 2 x 10 4 K/cm 3 , triangles 10 5 K/cm 3 , and squares 10 6 K/cm 3 . 
They are labelled according to the extinction Ay through the filament center. The dotted line in the left-hand figure 
shows the asymptote Eq. [17] for non-embedded cores {Ay = mag) , where the external radiation heating the core is 
completely absorbed. 



As seen in Fig. 1111 the curve provides an upper limit to 
the true L/M-ratios. For non-embedded cores pressurized 
by the mean ISM pressure, the ratio is, despite the high ex- 
tinction through the core center, lower by more then a fac- 
tor of two, even if the total luminosity of the SED (as given 
in Table [JJ) is considered. One reason is the steep density 
profile of the cores and the consequently strong variation 
in the optical thickness with the impact parameter of indi- 
vidual rays passing the core. A core that appears opaque 
in the center can still be optically thin or marginally opti- 
cally thick at the outskirts, so that only radiation towards 
the core center will be strongly attenuated. The amount of 
light absorbed in the core is further reduced by the clouds' 
albedo. 



The overall increase in the extinction with pressure 
(Ay oc y/p?) leads together with the decrease in the albedo 
to a stronger absorption in higher pressure regions. The 
L/M-ratio asymptotically approaches, as shown in Fig. [T2l 
the curve given by Eq. 1171 

As shown in Fig. [Til the simplified RT-model reproduces 
to first order not only the L/M-ratio but also the effective 
temperatures obtained using the full calculations. The de- 
viations are caused by the absence of the PAH emission 
in the simplified model, which is accurately accounted for 
in the full radiative transfer calculations. In the simplified 
model, all the absorbed UV and optical light is re-emitted 
in the FIR, which is otherwise re-emitted mostly by PAH 
molecules in the NIR and MIR. For extinction values below 
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Table 1. SED Parameters of Embedded Cores (T off = 10 K) 



A v 


Pc/k 




rji a 
dust 


/3 6 




L d 


L j A/ C orc 


T e 

-^total 


^total / -^corc 


FWHM f 


L Oauss 




[10 4 A] 

L cm J 


[%] 


[K] 




[—] 


[L Q ] 


'•AfpjJ 


[Le] 


L M. Q J 


[Rci] 


[-^inodel] 












(Pcxt = 


2 x 10 4 K/ 


™ 3 , X = 


1) 











2.00 


3.05 


18.05 


1.748 


2.697 


0.829 


0.323 


1.081 


0.421 


1.081 


1.077 










Cores in the center of 


spherical clouds 








1 


2.48 


2.92 


17.37 


1.737 


2.561 


0.562 


0.244 


0.675 


0.293 


1.059 


1.062 


2 


3.82 


2.79 


16.67 


1.714 


2.465 


0.343 


0.184 


0.391 


0.210 


1.062 


1.064 


4 


8.65 


2.83 


15.56 


1.682 


2.299 


0.142 


0.115 


0.152 


0.123 


1.086 


1.073 


8 


26.36 


3.02 


13.99 


1.627 


2.084 


3.99(-2) 


5.61(-2) 


4.09(-2) 


5.79(-2) 


1.131 


1.085 


16 


93.22 


3.51 


12.07 


1.549 


1.845 


8.12(-3) 


2.16(-2) 


8.13(-3) 


2.18(-2) 


1.176 


1.090 


DO 

32 


351. 1 


3.08 


y. /8 


1.553 


1.946 


1.37(-3) 


7.07(-3) 


1.39(-3) 


1 .19^-3 J 


1.22o 


1.101 


64 


1365. 


2.51 


7.66 


1.693 


2.637 


2. 30 (-4) 


2.34(-3) 


2.35(-4) 


2 39(-3) 


1.247 


1.105 










Cores in the center of c 


ylindrical clouds 








1 


2.72 


2.84 


17.09 


1.726 


2.561 


0.489 


0.222 


0.577 


0.262 


1.056 


1.062 


2 


4.76 


2.78 


16.26 


1.707 


2.443 


0.263 


0.158 


0.293 


0.176 


1.066 


1.066 


4 


12.46 


2.88 


14.98 


1.667 


2.256 


9.30(-2) 


9.04(-2) 


9.30(-2) 


9.57(-2) 


1.099 


1.076 


8 


42.04 


3.18 


13.28 


1.604 


2.025 


2.27(-2) 


4.05(-2) 


2.32(-2) 


4.14(-2) 


1.151 


1.090 


16 


158.1 


3.79 


11.30 


1.523 


1.801 


4.22(-2) 


1.46(-2) 


4.27(-3) 


1.48(-2) 


1.202 


1.098 












(Pcxt = 


1 x 10 b K/ 


cm 3 , x = 


1) 











10.0 


3.06 


17.28 


1.691 


2.329 


0.243 


0.211 


0.3021 


0.263 


1.151 


1.091 


Cores in the center of spherical clouds 


1 


10.48 


2.96 


16.68 


1.683 


2.298 


0.191 


0.170 


0.223 


0.199 


1.130 


1.083 


2 


11.92 


2.92 


16.17 


1.673 


2.250 


0.147 


0.139 


0.164 


0.156 


1.123 


1.082 


4 


17.33 


2.92 


15.26 


1.654 


2.161 


8.34(-2) 


9.56(-2) 


8.90(-2) 


0.102 


1.128 


1.084 


8 


37.01 


3.10 


13.84 


1.608 


1.996 


3.02(-2) 


5.05(-2) 


3.11(-2) 


5.20(-2) 


1.154 


1.090 


16 


108.7 


3.56 


12.08 


1.530 


1.767 


7.06(-3) 


2.03(-3) 


7.14(-3) 


2.05(-2) 


1.202 


1.104 


32 


378.2 


3.64 


9.78 


1.528 


1.838 


1.25(-3) 


6.70(-3) 


1.27(-3) 


6.81(-3) 


1.233 


1.103 


64 


1417. 


2.57 


7.64 


1.662 


2.463 


2.10(-4) 


2.18(-3) 


2.14(-4) 


2.22(-3) 


1.262 


1.110 










Cores in the center of c 


ylindrical clouds 








1 


10.73 


3.00 


16.49 


1.678 


2.315 


0.178 


0.160 


0.203 


0.184 


1.122 


1.081 


2 


12.89 


2.96 


15.86 


1.666 


2.251 


0.126 


0.124 


0.138 


0.137 


1.118 


1.081 


4 


21.14 


2.97 


14.75 


1.644 


2.159 


6.22(-2) 


7.87(-2) 


6.55(-2) 


8.30(-2) 


1.131 


1.085 


8 


52.21 


3.24 


13.20 


1.587 


1.954 


1.90(-2) 


3.78(-2) 


1.94(-2) 


3.86(-2) 


1.161 


1.088 


16 


171.9 


3.82 


11.25 


1.521 


1.784 


3.92(-3) 


1.41 (-3) 


3.97(-3) 


1.42(-2) 


1.208 


1.099 












(Pcxt = 


1 x 10 fa K/ 


™ 3 , X = 


1) 











100.0 


3.79 


16.48 


1.498 


1.480 


3.47(-2) 


9.57(-2) 


0.0409 


0.113 


1.311 


1.125 


Cores in the center of spherical clouds 


1 


100.5 


3.70 


15.87 


1.505 


1.541 


2.94(-2) 


8.13(-2) 


3.29(-2) 


9.09(-2) 


1.283 


1.111 


2 


102.0 


3.59 


15.40 


1.510 


1.571 


2.53(-2) 


7.02(-2) 


2.74(-2) 


7.62(-2) 


1.280 


1.119 


4 


107.8 


3.47 


14.62 


1.516 


1.605 


1.89(-2) 


5.39(-2) 


1.98(-2) 


5.67(-2) 


1.256 


1.108 


8 


130.4 


3.50 


13.43 


1.508 


1.607 


1.08(-2) 


3.38(-2) 


1.10(-2) 


3.47(-2) 


1.245 


1.105 


16 


214.4 


3.83 


11.82 


1.467 


1.538 


3.97(-3) 


1.60(-2) 


3.99(-3) 


1.62(-2) 


1.247 


1.106 


32 


518.7 


3.54 


9.644 


1.501 


1.709 


9. 26 (-4) 


5.81(-3) 


9.42(-4) 


5.91(-3) 


1.278 


1.121 


64 


1638. 


2.68 


7.60 


1.607 


2.181 


1.71 (-4) 


1.91(-3) 


1.75 (-4) 


1.95(-3) 


1.298 


1.123 










Cores in the center of c 


ylindrical clouds 








1 


100.7 


3.68 


15.58 


1.518 


1.637 


2.83(-2) 


7.83(-2) 


3.13(-2) 


8.66(-2) 


1.279 


1.121 


2 


102.9 


3.59 


15.04 


1.520 


1.654 


2.33(-2) 


6.51(-2) 


2.50(-2) 


6.98(-2) 


1.254 


1.107 


4 


111.7 


3.51 


14.16 


1.522 


1.674 


1.62(-2) 


4.72(-2) 


1.69(-2) 


4.92(-2) 


1.241 


1.104 


8 


145.8 


3.61 


12.87 


1.504 


1.643 


8.22(-3) 


2.73(-2) 


8.38(-3) 


2.79(-2) 


1.238 


1.104 


16 


275.1 


4.01 


11.12 


1.473 


1.609 


2.62(-3) 


1.20(-2) 


2.65(-3) 


1.20(-2) 


1.246 


1.105 



Effective dust temperature of the modified black-body fit. 

Power of the effective emissivity kT = «:o(A/A )- /3 where Ko is the effective emissivity at Ao = 250/^m. 



Effective emissivity per dust mass at Ao = 250 (jm. 
Luminosity of the modified black-body fit. 
Total luminosity of the dust re-emission. 

FWHM of the Gaussian source approximation of the surface brightness at 250 fim. 
Flux of the Gaussian source approximation of the surface brightness at 250 /im. 



Ay = 4 mag, the heating of the dust grains in the cores is 
overestimated leading to a higher dust temperature. In the 
case of highly opaque filaments (above Ay ~ 8 mag), the 
PAH emission of the filaments leads in the full RT calcu- 



lation to an increasingly important heating source, which 
gives rise to the warmer dust temperatures (Sect. [A"l). 

The curves for the lowest external pressure provide a 
lower limit to both the L/M-ratio and the dust tempera- 
ture of cores. Above the curves, passively heated cores are 
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Fig. 10. Relation of the effective values kq and /3 derived 
by fitting a modified black-body function to the theoreti- 
cal dust emission spectrum of condensed cores. The sym- 
bols represent the same external pressure p ex t/k of the 
medium surrounding the filaments as in Fig. [S] (diamonds 
for 2 x 10 4 K/cm 3 , triangles for 10 5 K/cm 3 , and squares 
for 10 6 K/cm 3 .) The symbols corresponding to filaments 
pressurized by the mean ISM pressure are labeled with the 
assumed central extinction of the filament. The curves are 
the derived fits for (3 and kq shown in Fig. [9j 
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Fig. 12. Light-to-mass ratio (L/M) of the dust luminosity 
and the core mass as function of the dust temperature of 
embedded cores. The derived ratios are labelled with the 
assumed extinction values Ay through the filaments cen- 
ter. The ratios L / M corc (T<iust) obtained for cores that are 
embedded in cylindrical and spherical filaments for each 
given external pressure p ay ±jk (2 x 10 4 , 10 5 , 10 6 K/cm 3 ) 
are fitted by a polynomial function shown as solid black 
lines. The points for mean ISM pressure are fitted also by 
a simple power law (dotted curve). 



embedded in filaments, which are surrounded by a higher 
pressure gas. The cores are also located in this region if they 
are heated by an additional source such as a progenitor star 
in the core or an additional external source. 

The resulting relation between the L/M-ratio and the 
effective dust temperature Tdust for a given external pres- 
sure is shown in Fig.[TJ] For the same L/M-ratio, the dust 
in embedded cores appears to be warmer in high pressure 
regions of the ISM because the SED is largely determined 
by the relative warm dust emission at the outskirts of the 
cores. The L/M-ratios for low external pressures lie on a 
curve that is slightly flatter than a simple power law as 
shown in the figure. The power-law fit provides a constant 
emissivity with (3 — 1.96 and kq — 2.15 cm 2 /g. 

3.3. Brightness profile 

Radiative transfer effects cause the brightness profiles of 
the dust emission spectrum to be broa der than the intrin- 
sic pr ofile of the column density profile (jFischera fc Dopital 
2008). The profiles of the surface brightness at three far 
infrared wavelengths, which correspond to the bolometers 
of the Herschel satellite's SPIRE-instrument (and of the 
balloon experiment BLAST), are shown in Fig. [T3J The 
profile steepens towards longer wavelengths. This effect is 
caused by the emission at longer wavelengths being dom- 
inated by larger grains, which have lower radial tempera- 
ture variations than the smaller dust grains that dominate 
at shorter wavelengths. In addition, the emission at shorter 
wavelengths is more affected by temperature variations be- 
cause these regions are blueward of the Rayleigh- Jeans 
limit. The steepening of the profile increases for highly em- 
bedded cores (as the cores are more optically thick) that 
are affected by stronger radial temperature variations. 



The cores are often modeled obs e rvatio nally as Gaussian 
source s (see e.g., iNetterfield et al.l (|2009f ) and lOlmi et ail 
(2009])). To facilitate comparison of observations to my 
models, I approximate the surface brightness at 250 /^m 
by a two-dimensional Gaussian function. A Gaussian func- 
tion is broader than the profile of the column density, but 
as can be seen in Fig. 1131 a Gaussian function can be a 
good representation of the true emission profile as in case 
of cores embedded in filaments with Ay < 8 mag that are 
confined by the mean ISM pressure. The FWHM of the 
Gaussian function is a good estimate of the physical radius 
of the core (Fig. [TJ} albeit too large by 10%. The Gaussian 
approximation becomes worse in high pressure regions or 
for highly embedded cores where the cores are more com- 
pact. The error in the size estimate based on the FWHM 
increases therefore towards lower dust temperatures. The 
Gaussian approximation for highly embedded cores is too 
narrow, overestimating the flux passing through the core 
center and inferring a FWHM considerably larger than the 
intrinsic core radius (Fig. RH)) . The estimated flux using a 
Gaussian approximation is slightly too high (6% to 12% 
overestimate) depending on the external and internal pres- 
sure of the filament. The flux estimate tends to be more 
accurate in the lower pressure regions of the ISM. 

4. Discussion 

My model connects the effective dust temperature to the 
physical properties of passively heated embedded cores. In 
particular, it establishes a relation between the dust tem- 
perature and the pressure of the medium surrounding the 
cores. The details of the relationship are determined by the 
ISRF, the dust properties in the cores, and the external 
pressure confining the filaments but is independent of the 
actual mass of the cores. This model can be directly tested 
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Table 2. Polynomial approximations 



[K/cm 3 ] 


do 


«i 


02 


a3 


tt4 








P 


= Ei a; (log 


io TdustfK]) 1 






2 x 10 4 


-40.7509 


186.611 


-296.453 


202.709 


-50.5764 


— 


1 x 10 5 


-41.0076 


190.600 


-307.465 


213.560 


-54.1691 


— 


1 x 10 6 


-29.0751 


144.608 


-243.478 


175.672 


-46.2471 


— 






kq [cm 


7g] = ^ 


(log 10 T dust [K])* 






2 x 10 4 


-133.777 


600.783 


-948.039 


639.586 


-156.646 


— 


1 x 10 5 


-142.593 


654.894 


-1060.36 


735.858 


-185.995 


— 


1 x 10 6 


-90.8491 


461.490 


-802.533 


592.289 


-158.761 


— 






log 10 L/M corc [L Q /M Q ] 


= Eidi(log 10 Tdust 


[K]y 




2 x 10 4 


-2.63348 


-4.06995 


4.60187 








1 x 10 5 


-3.27281 


-2.87834 


4.03329 








1 x 10 6 


-6.11354 


2.87895 


1.09179 












FWHM[R cl ] = 


E iai (T dust [K]) 1 






2 x 10 4 


1.04042 


0.193172 


-0.0506557 


0.00579236 


-0.000312315 


6.33987e-06 


1 x 10 5 


1.01551 


0.198970 


-0.0502195 


0.00571245 


-0.000313287 


6.59132e-06 


1 x 10 6 


1.29788 


-0.0283919 


0.0144707 


-0.00232888 


0.000140731 


-2.79548e-06 
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Fig. 14. FWHM at A = 250 ^m of a Gaussian represen- 
tation of the surface brightness profile of cores which are 
either embedded in spherical or cylindrical filaments. The 
data points for each assumed external pressure are fitted 
with a polynomial function shown as solid lines. 

by high resolution observations in the FIR/submm regime, 
which allow apart from measuring the core mass by model- 
ing the dust emission spectrum also a determination of the 
size of the core and by applying the mass-radius relation 
(Eq. [3]) the pressure around the core. 

The model relates essential core parameters such as the 
FWHM, the emissivity (/3 and kq), and the luminosity- 
mass ratio to the pressure of the medium surrounding the 
filaments and the extinction through the filaments' center, 
which corresponds (for a given filament shape) to a certain 
central pressure. The values of the parameters are found 
to be (for a given external pressure) controlled mainly by 
the dust temperature, which allows in principle a simple 
correction for radiative transfer effects. Future work needs 
to determine how shapes of the ISRF, the dust properties, 
and the structure of the molecular clouds caused by the 
turbulent motion affect the core models, allowing a clearer 
understanding of the uncertainties associated with the mea- 



surement of core masses. The expected basic effects are dis- 
cussed in Sects. I4.2.H 14.2.21 and 14.31 A additional compli- 
cation is caused by the coupling of gas and dust in high 
density environment, although the effect might be less im- 
portant if the cores are mainly supported by turbulent mo- 
tion. 

4.1. The origin of the effective dust emissivity 

The emissivities of my models differ considerably from the 
mean properties derived in Sect. l2.2?T1 There are two major 
reasons why the dust emissivity seems to be less steep and 
also lower. First, the variation in dust temperatures with 
size and composition at each radial distance in the core and 
second the smooth variation in the dust temperatures with 
position in the core. 

Even without additional radiative transfer effects the 
parameters of the emissivity have been shown (Sect. I3.2[) 
deviate from the intrinsic mean values. The reason is that 
the emission behavior depends on both grain size and grain 
composition. The mass of the dust grains is contained in the 
large grains, which are also the coldest with temperatures 
much below 19 K. The warm dust temperature is caused by 
the smallest dust grains, which do not contain much mass 
but efficiently absorb interstellar radiation. As a result their 
emission dominates the blue side of the SED and, to a large 
extend, its peak. At larger wavelengths, the overall dust 
emission is gradually enhanced by the emission of colder 
and larger grains, which flattens the emission profile. To 
compensate for the warm dust and the low mass related to 
the dust emission, the emissivity is reduced to obtain the 
true dust mass. 

A similar effect is caused by the gradually changing dust 
temperature with radius. The emission from warm dust 
grains located at large cloud radii dominate the bulk of the 
SED and determine to a large part the position of the peak 
of the dust emission and therefore the effective dust tem- 
perature. This is most clearly demonstrated by the similar 
effective temperature shown in Fig. [11] and the dust tem- 
perature at the cores' edge as shown in Fig. [5J At longer 
wavelengths, the integrated dust emission spectrum is grad- 
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Fig. 13. Brightness profiles at 250, 350, and 500 /im. They are compared with the profile of the column density (dotted 
line). The surface brightness at 250 /jm is fitted to a Gaussian source (grey line). The filaments arc cither confined by 
a medium with p ex t/k = 2 x 10 4 K/cm 3 (top) or 10 6 K/cm 3 (bottom). The extinction through the filaments' center is 
either Ay — 8 mag or Ay — 64 mag (left-hand figure and right-hand figure, respectively). 



ually more affected by colder dust grains in the interior of 
the cloud. The shape of the dust re-emission overall appears 
to be flatter than expected for the underlying dust emissiv- 
ity As the dust emission is dominated by the warm dust 
grains, the absolute emissivity needs to be reduced to ac- 
count also for the emission of cold dust grains. As seen, the 
effect can be larger than a factor of two depending on the 
radiation field heating the core and how opaque the core 
actually is. 

As seen, the effect of the temperature distribution inside 
the core on the apparent emissivity increases for a given 
central extinction of the filaments with external pressure 
and (for given external pressure p c *t/k < 10 6 K/cm 3 ) with 
extinction of the filaments until the filaments become highly 
opaque with Ay > 16 mag. The steepening of the apparent 
emissivity for embedded cores in filaments with Ay > 16 is 
caused by an overall reduction in the temperature variation 
inside the core both in terms of the location as seen in 
Fig. [5] and in terms of grain composition and grain size that 
occurs for of highly reddened radiation fields. The reason for 
the latter phenomenon is given in Sect. 13.1.31 The reduced 
radial temperature variation as seen in Fig. [3] is caused by 
the additional heating of the grains in the core center by the 



dust emission from the filament. A similar effect is expected 
if the dust emission from the Galactic disc is considered. 

A physical limit exists for embedded cores that become 
predominantly heated by the CMB background radiation. 
In this situation, the dust in the core approaches the tem- 
perature of the CMB independently of the position and 
composition or size and the apparent emissivity becomes 
equal to the intrinsic dust emissivity. 

The temperature variation is probably also the explana- 
tion of the flatter emissivity observed in other astrophysical 
objects such as galaxies for e xample (|Dunne fc Ealesll200TI: 
IChakrabarti fc McKeeH2008l) . However, in the literature the 
possibility of a dependence of the emissiv ity on the dust 
prope rties is explored (see discussion by iDunne fc Ealesl 

mm). 

4.2. Basic model uncertainties 
4.2.1. Uncertainties in the ISRF 

The ISRF applied in the model is a canonical approach to 
calculate the dust re-emission from dust grains located in 
the ISM. However, a large fraction of the radiation in the 
galactic disc is caused by dust grains and PAH molecules, 
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Fig. 15. Mean interstellar radiation field assumed to 
model the SEDs of embedded cores. The solid black line 
is the ISRF c aused by the stellar light as provided by 
iMathis et al.l (|l982f) . The grey line shows the assumed 
angle-averaged dust re-emission. The contribution of the 
PAH emission is reduced relative to the dust re-emission 
of the diffuse ISM (Fischera, in prep.) shown as light grey 
curve. The strength of the radiation field is scaled to the 
DIRBE flux at 140 ftm. The DIRBE fluxes at 100, 140, and 
240 /im are shown as open diamonds. The dust emission 
is divided into its several components: PAH (long dashed 
line), graphite (short dashed line), silicates (dotted line), 
and iron (dashed-dotted line). 



which efficiently absorb the light in the UV and optical 
light. For comparison with the results presented in this pa- 
per, I considered an ISRF where the dust re-emission from 
the Galaxy is taken into account. The most important con- 
tribution arises from the PAH emission and possibly small 
stochastically heated grains. In general, this problem needs 
to be addressed more quantitatively based on all sky sur- 
veys such as the DIRBE experiment. The determination 
of the corresponding mean intensity inside the Galaxy is 
complicated by the spectral range being dominated by the 
zodiacal light. The approach presented here is qualitative, 
but demonstrates where the heating by dust grains and 
PAH molecules becomes important. 

The mean intensity of the dust emission consists pri- 
marily of three components: emission from the diffuse gas, 
HH-regions, and dense clouds (which potentially contain 
a large percentage of the interstellar material consider- 
ing that prob ably up to 30% of hydrogen is molecular 
(|FerriereJ 12001 ). One of the main differences between the 
three components is the strength of the PAH emission. In 
dense clouds, it is likely that the PAH molecules are con- 
densed onto grains. If present, the PAH molecules would 
be heated by highly attenuated optical and UV radiation, 
which would reduce considerably the PAH emission rela- 
tive to the dust emission peak. Additional PAH emission 
and warm dust emission (around 60 K) arises from HII- 
regions. I assume that the average PAH emission is lower 
than the PAH emission from the diffuse phase by a fac- 
tor of 2. Otherwise the dust re-emission is assumed to be 
the same as in the diffuse ISM. I use the same dust prop- 
erties as in the ray-tracing calculations. The SED of the 
dust emission is derived assuming a non-attenuated ISRF. 
The strength of the dust emission is obtained by scaling 



the emission at 140 /im to the corresponding mean inten- 
sity of the DIRBE experiment (yj v = 6.78 x 10~ 7 W/m 2 ). 
The assumed total SED of the mean ISRF is displayed in 
Fig. HS1 I derived the DIRBE fluxes shown in the figure us- 
ing all sky DIRBE data from which the zodiacal light had 
been subtracted using a physical model. The colder dust 
emission from compact clouds relative to the emission from 
diffuse dust grains has only a small effect on the heating 
of the cores and is ignored. I also neglect the warm dust 
component from HII- regions. In this model (Fig. IT51) . 34% 
of the Galactic mean intensity in the solar neighborhood is 
caused by dust grains and PAH molecules. 
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Fig. 16. Effective temperatures of embedded cores as a 
function of the central pressure in the filament, which is 
assumed to be either a self-gravitating sphere or a self- 
gravitating cylinder. In case a), the emissivity is assumed 
to be the same as in the diffuse dust model, while in 6) the 
emissivity is increased by a factor 8. The external pressure 
outside of the filaments is assumed to be 2 x 10 4 , 10 5 , and 
10 6 K/cm 3 . The filaments are h eated either by the ISRF as 
provided bv lMathis et al.l (|l982l ) (grey lines) or by an ISRF 
where the diffuse dust emission is added. 

In Fig. [T51 the PAH emission of the assumed ISRF is 
clearly stronger (despite it being fainter than the diffuse 
dust emission) than the PAH emission at the center of the 
filaments (Fig. [3]). This implies that, even if PAH emission 
arises from the filaments, it is less important for heating 
the dust in the core than the Galactic IR emission. The 
result for the modified ISRF is shown in Fig. [T5] The ad- 
ditional heating by IR photons becomes important only for 
highly embedded cores where the central extinction of the 
filaments is larger than Ay ~ 8 mag. This additional heat- 
ing is more important for the inner core regions, which lead 
to a steeper brightness profile. 

4.2.2. Uncertainties in the emissivity 

The results presented in this paper depend to a certain de- 
gree on the dust model, which determines both the position 
of the peak of the dust rc-cmission spectrum and its shape. 
A different dust model will therefore lead to different ef- 
fective dust parameters being assigned to the derived dust 
emission spectra for embedded cores. This will also be true 
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if the results a re compared with calcu l ations based on the 
dust model of IWeingartner fc Draind (|200lf) . which has a 
different assumed dust composition and size distribution. 
However, since both models were constructed to reproduce 
the observed properties in the diffuse ISM the differences 
should be small relative to the large uncertainties in the 
dust properties of the dense phase. 

The calculations are based on the dust properties of the 
diffuse medium merely as a reference point. Grains in the 
dense phase are likely to form icy layers that will, as a 
consequence, increase the sticking coefficient and therefore 
accelerate the grain growt h process by means o f grai n coag- 
ulation (see discussion in iFischera fc Dopital ()2008[) V The 
process will naturally alter the size distribution of the grains 
by removing the number of small grains while increasing 
the mass of the largest dust particles. This will result in a 
rather flat slope of the extinction curve in the UV, which 
is generally determined by the abundance of small grains. 
Its effect on the optical properties in the far infrared is 
uncertain but the emissivity may be c o nside rably higher. 
According to iKruegel fc Siebenmorgenl (|1994[ ) , a factor of 
10 higher emissivity is physically possible depending on the 
properties of the co mposites. In addition, the grain growth 
models presented by lOssenkopf fc Henningj (|1994l ) depend 
on the density inside the clouds. If correct, a higher emis- 
sivity is expected in denser clouds. The mass estimate of 
cores is observati onally based on an emissivity at 250 /xm 
of k = 10 cm 2 /g (jNetterfield et al.ll2009l: lOlmi et al.ll2009D 
with (3 = 2 and a gas-to-dust mass-ratio of 100. Despite 
this relatively low gas-to-dust mass ratio, the cores appear 
to be far less massive than they would be if the masses were 
based on the emissivities derived in this paper. 

Figure shows a first-order effect of a higher emis- 
sivity on the core temperature. It is assumed that the 
spectral index j3 has the same dependence on the pres- 

(2) 

sure. In this case, the core temperature varies as T^ u ' st = 

-^dust( K o,i/ K o,2) 1 ^ 4+ ' 3 ^ In the special case where the emis- 
sivity of the dust grains has increased by a factor 8, the 
absorption coefficient above 60 /im is scaled by the same 
factor. Between 20 /im and 60 /Ltm, a smooth transition is 
made. The pressure p c in the filaments obtained by mea- 
suring the mass by modeling the SED and the size of the 

core varies as = Pc (^o,i/ K o,2) 2 and therefore depends 
strongly on the emissivity. Too low emissivities will produce 
too high pressure values. The uncertainty in the emissiv- 
ity affects observationally only the mass estimate and the 
pressure. The P — T-relation might therefore be helpful to 
estimate the dust emissivity in cores. 

A first-order correction to the dependence of /3, kq, 
L/M corc , and FWHM on Tdust is given by replacing the 

(2) 

dust temperature with Xj'*. For the emissivities kq, the 
values have to be corrected by the corresponding factor. 
However, the dependence on dust temperature needs to be 
analyzed quantitatively using more realistic assumptions 
about the dust properties inside the clouds. 



4.3. Effect caused by the inhomogeneous medium 

As mentioned in the introduction, molecular clouds are tur- 
bulent, which produces a strong variation in the local den- 
sity and as a consequence a variation in t he column density 
seen t hrough the molecular clouds (e.g. Kauffm ann et al.1 
(|2010h ). The effect of turbulence was studied in some de- 



tai l bv IFischera fc Dopital (12004 l2005f ) and discussed also 
bv IFischera fc Dopita (|2008T ). where a method is described 
that considers how the turbulent density structure in ra- 
diative transfer studie s through self-gravitating clouds can 
be taken into account. [Fischera fc Dopital (|2005l ) illustrated 
how the dust attenuation can be directly related to the 
physical condition of the turbulent medium. A statistical 
treatment of the radiative transfer through the clouds show- 
ing a smoothly varying density fluctu ations as expected 
in tur bulent clouds was presented by fHcgman n fc Kegel 

4111). 

Despite a number of studies about the effect of 
an inh o moge n eous medium o n the radiative t r ansfer 
(iBoissd 119901: IWitt .fc Gordon! 119961: IWolf et all l l998: 



I Varosi fc Dweklll999l: IWitt fc Gordo"nTl2000 | : IGordon et all 
l200lt iHegmann fc Kegel! 120031: iDotv et al.ll2005l ). the im- 

portance for the interpretation of infrared observations is 
still uncertain. Although similar in nature, not all stud- 
ies can be directly applied to the problem of externally 
heated molecular clouds. The problem generally addressed 
in studies about the radiative transfer through molecu- 
lar clouds is the effect of the density structure on the 
dust scattering and d ust absorption without considerin 



the dust re-emission dBoissd [l990t IWitt fc Gordon! Il99' 



Wolf et al.lll99l IWitt fc Gordon! ^OOOU Hcgman n fc Kegel 
120031) . While it is one of the main complica tions of radia- 
tive transfer calculations, IDotv et ail (|2005l ) stated (based 
on test calculations) that dust scattering plays only a very 
small role in the distributions of dust temperatures in the 
majority of their sources. They indeed ignored this problem 
altogether in their study of dust temperatures in molecular 
clouds. 

Because of the complexity of the problem, the model 
calculations have in most cases been based on rather sim- 
plified assumptions about the medium. In general, a two- 
phase model is considered where randomly distributed 
dense 'clumps' are e mbedded in an otherwise homoge- 
ncous thin medium dBoissd 119901: IWitt fc Gordon! 119961: 
Wol f et al.l[l998t IWitt &: Gordonll2000l: IGordon eHal boOl: 
IDotv et aLll2005l ). The inhomogeneous medium is typically 
described by the filling factor of the 'clumps', their size rel- 
ative to the extension of the whole cloud, the density ratio 
of the 'clumps' and the 'inter-clump' medium, and in ad- 
dition by the optical depth of the cloud. The application 
to observations is complicated by neither the medium as a 
whole nor the individual clumps being described as physical 
entities with realistic density profiles. 

To some extent, the models of individual dense 'clumps' 
are comparable to the dense self-gravitating cores in my 
model of embedded cores. From a physical standpoint, it 
seems, however, more realistic that they form in the inner 
center of the filaments where they experience the highest 
pressure rather than to place them randomly inside the 
filament. Because of the pressure profile, it is also expected 
that cores located closer to the cloud edges are less optically 
thick than at the center. My model is less applicable to 
cores formed in a dense cluster where shadow effects by 
individual cluster members become important. Although 
this situation is often considered in two phase models by 
applying a high filling factor of the 'clumps', it becomes 
observationally difficult to separate the individual sources. 

The most important characteristic caused by in- 
homogeneities is that a non-homogenous dusty medium 
is more transparent than an inhomogeneous medium 
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(iBoissd I1990L iHegmann fc Kegel! l2003t iDotv et al.l 120051: 
iFischera et all 120031: iFischera fc Dopital I2005T I . The effect 
disappears when the density enhancements become opti- 
cally thin. For typical molecular clouds (Ay ~ 8 mag), the 
increase in transparency is highest in the UV and optical 
and is negligible in the infrared. The effect leads there- 
fore for a homogeneous extended source seen through a 
turbulent screen to a fla ttening of the attenuation curve 
(|Fischera fc D opita 2005). This flattening is characterized 
by a large absolute-to-relative extinction, the i?y-value be- 
ing defined as Ay /(Ab — Ay), which is close to 3.1 in the 
diffuse ISM. The same effe ct is expected for th e turbu- 
lent interstellar medium and IFischera et al.l (|2003h demon- 
strated that the turbulence can naturally explain the flatter 
curvature of the Calzetti extinction curve derived for local 
starburst galaxies. Because of the inhomogeneities, the light 
most responsible for heating the gas and the dust can there- 
fore penetrate deeper into the filaments and consequently 
into the core. On the other hand, the total light absorbed 
inside the core is lower leading effectively to a lower dust 
temperature, solely by considering the energy balance in 
the core. 

In the framework of the model, the SEDs of the cores are 
independent of the effective temperature, which determines 
for a given pressure the mass and the size of the cores in- 
side the filaments. If the cores are supported predominantly 
by turbulent pressure, more massive cores are also more 
turbulent, which creates a stronger variation in the local 
density around the mean. Their one-point statistic of the 
local density is described by a broader probability distri- 
bution approximately described by a log-normal function. 
It is therefore expected that in a certain pressure region 
of a GMC, the effect of the turbulent density structure is 
stronger for more massive cores than for less massive ones. 

5. Summary and conclusions 

I have presented a physical model for the dust emission 
of condensed cores where the cores are described by crit- 
ical stable self-gravitating clouds embedded at the center 
of self-gravitating filaments. The filaments are assumed to 
be either spherical or cylindrical in shape. For a set of dust 
properties and an external isotropic radiation field, the ra- 
diative transfer problem is determined by the external pres- 
sure confining the filaments and the extinction of the fila- 
ments. In the model, the extinction of either a sphere or a 
cylinder is directly related to the overpressure in the fila- 
ment center and the external pressure. 

The problem has been solved using a ray-tracing tech- 
nique where the effects of multiple scattering events and 
the anisotropy of the scattered light are accurately taken 
into account. The dust emission is derived using a phys- 
ical dust model of stochastically heated grains of differ- 
ent sizes and compositions. The work focuses on the ra- 
diative transfer effects and is therefore based on the dust 
properties of diffuse dust grains. The mean extinction at 
FIR wavelengths is approximately described by a simple 
power law with /3 cxt = 2.024 and an absolute extinction at 
A = 250 fim of K oxt (A ) = 5.232 cm 2 /g- 

The variation in the dust temperatures inside embed- 
ded cores is found to depend on the external pressure and 
the extinction of the filaments. The variation is not nec- 
essarily smaller than for non-em bedded cores as stated by 
IStamatellos fc Whitworthl (|2003h . 



For comparison with observations, the bulk of the dust 
re-emission is approximated by a modified black-body spec- 
trum using /3 and kq of the emissivity k™ 1 = Ko(A/Ao) - ^ f° r 
Ao = 250 /im and the dust temperature Tdust as free param- 
eters. As shown, the effective parameters of the emissivity 
can be applied using a single-zone model for the cores. This 
simplified model reproduces, to first order, the dust temper- 
ature of the cores. On the basis of the model assumptions 
of self-gravitating clouds, the dust temperature of passively 
heated cores is independent of the mass of the cores. It is 
found that: 

— The effective emissivity kq is considerably lower than 
the intrinsic extinction. For filament extinction Ay < 
16 mag, the values of kq are typically a factor 2 to 3 
lower than expected using the mean dust properties in 
the cloud. 

— The emissivity is flatter with (3 < 1.8. For the mean 
value of ISM pressure and Ay < 16, 1 found that 1.52 < 
P < 1.74. 

— The effect on kq and /3 for a given dust temperature 
increases with external pressure. 

— In highly embedded cores with Ay > 16, the emissivity 
becomes equal to the mean extinction value. 

The derived surface brightness is considerably broader 
than the underlying profile of column density, the effect 
being most prominent at shorter wavelengths. The bright- 
ness profiles become flatter for both more embedded cores 
and higher pressures of the external medium. The surface 
brightness at 250 /im has been compared with a Gaussian 
source approximation. The agreement is good for less em- 
bedded cores, but is poor for strongly embedded cores, 
where the Gaussian profile is steeper than the model. It 
is found that: 

— The Gaussian source approximation overestimates the 
total flux by ~ 10%. 

— In high pressure regions and highly embedded cores, the 
FWHM of the Gaussian source approximation overes- 
timates the size of the core by between 10% and 30%. 

For direct comparison with observations I have derived, 
for different ISM pressures, polynomial expressions of the 
model relations for /3(T dust ), K (T dust ), V-^core^dust), and 
FWHM(T dust ) at A = 250 /jm. The relations are found to 
be only mildly affected by the geometry of the filaments. 
For highly embedded cores, the work indicates that the dust 
emission from the Galactic disc needs to be considered to 
model accurately the SED and surface brightness profiles 
of the cores. 
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Appendix A: Approximate solution for the dust 
temperature in cores 

In a single temperature model used by observers to inter- 
pret their observations and in this paper, the individual 
dust grains are, independently of their location inside the 
core, heated by the same mean intensity J^ cat . In the frame- 
work of the model, this intensity is produced by the radia- 
tion field at the center of the filament and by the emission 
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from the dust grains in the core. That the dust temperature 
does not strongly depend on the heating flux (it varies as 
the sixth root of the heating rate) allows approximations 
that I outline in the following. The aim of the model is 
to provide an equivalent to the time-consuming ray-tracing 
technique. To be widely applicable, the effects caused by 
scattering and re-heating by dust grains are taken into con- 
sideration. As in case of the more accurate approach using 
ray-tracing, the solution is derived in two independent con- 
secutive steps by estimating first the mean intensity J\(0) 
in the center of the filament and secondly the dust temper- 
ature inside the core. 

To obtain a most accurate estimate of the attenuated 
flux caused by the self-gravitating cores and filaments hence 
a more reliable estimate of the absorbed energy inside the 
clouds, I still consider clouds with physical density profiles 
rather than homogeneous clouds. The difference between 
the important optical properties (their attenuation and es- 
cape probability) for the two assumptions about th e den sity 
distribution inside the clouds is discussed in Sect. IB .11 and 
IB. 21 where it is shown that an assumption of a homoge- 
neous cloud will lead to an over-estimate of the absorbed 
flux. 

In both steps, the dust temperature in the correspond- 
ing cloud is derived. To consider the radiative transfer ef- 
fects caused by the emission from dust grains, it is assumed 
that the emissivity r^ ust = hf* B x (T duat ) in both the fila- 
ment and the core is proportional to the hydrogen density 
n H . 

There are in general three different contributions to the 
intensity J\(0) at the center of the filament heating the 
core: (1) the attenuated radiation of the interstellar radia- 
tion field illuminating externally the filament, (2) the radi- 
ation of the scattered emission, and (3) the dust emission. 
To derive accurately the last two contributions, a proper 
radiative transfer treatment is needed. 

The light scattered away from the initial direction is 
partly compensated for by radiation from other initial di- 
rections scattered into it. This leads effectively to a lower 
extinction towards the cloud center. The approximation 
where the light is scattered in either a forward or back- 
ward direction and where the corresponding fraction is de- 
termined by the g-factor, the angle-averaged cosine of the 
scattered intensity, leads to a lower extinction 
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where uj\ — fc A ca /fc^ xt is the dust albedo. In the UV and 
optical, the light is strongly scattered in the forward direc- 
tion with g\ > 0.5 and ui\ > 0.5. This causes the extinction 
at these wavelengths to be considerably lower than nomi- 
nal values. The assumption for scattered light has overall 
less effect for highly embedded cores where the contribution 
of the scattered emission to the mean intensity inside the 
filament is weaker. 

The dust grains that provide most of the intensity at 
the center of an optical thick filament are located at the 
edge of the filament where the dust temperatures are the 
warmest. Further inside, the grains are heated by a strongly 
attenuated radiation field, which shifts the dust emission 
spectrum to longer wavelengths, where the absorption is 
also weaker. A conservative lower limit to the dust heating 
is provided by assuming that all the dust emission arises 



from a thin shell at the outer surface of the cloud. This 
provides 
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for a sphere and 

Jt ust (0) = \(vt ust /n^(N R )(0) (, - 
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for a cylinder, where (iVn) (0) is the corresponding column 
density of hydrogen at the center of the homogeneous cloud 
(Fig. IB. II) . The mean attenuation is given by 
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where the optical depth given by 
r'r' I ;/.-;r n H )n H (r)dr. 



(A.4) 
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The dust emission of the filament and the core are cal- 
culated using energy conservation between the absorbed 
and the emitted flux. The absorbed luminosity in the cloud 
caused by the external radiation is given by 



Lf s ' a = S d nJ x (l-exp{-Tf (p,r A )})(l-^ 1 (p,r A )) ) (A.6) 



where S c \ is the cloud surface, r^ ff the effective optical 
depth fsee lB.lJ) . and cof the albedo of the cloud f see IB. 3ft . 
The attenuation and the cloud albedo depend (apart from 
the shape) on the overpressure p — p c /p cx t and the ra- 
dial extinction t\. The additional dependence on the non- 
isotropy of the scattered light is ignored. The luminosity 
of the dust emission absorbed inside the cloud is obtained 
using the escape probability, p esc (see IB. 21) , 
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( 1 - Pose (P, T\ ) ) ( 1 - Ujf (p, T A ) ) , 
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where Lf ust = V c i (ran) 47r?7^ ust /^H, V c \ is the cloud vol- 
ume, and (jih) the mean hydrogen density. The absorbed 
luminosity is, on the other hand, equal to 



t abs 



V c i (n H ) (fcf7n H )47rJf at . 
The external radiation field is therefore reduced to 



J heat, a r / - \ T 
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with the mean extinction of the cloud defined as (t\) = 
4Vci/5 c i (tin) fc A bs /riH- The intensity produced by the dust 
grains emission is given by 



J 



hcat.b 



^dust 



/kf S (1 - Pesc (p, Tx))(l - LJ?(p, Tx)). (A.ll) 



The contribution of the thermal dust emission to the total 
heating rate from dust grains located in the cloud is derived 
in an iterative manner. 

The contribution of the different components of the ra- 
diation considered in the approximation to the total heating 
rate are shown in Fig. IA.1I for cores embedded in spher- 
ical filaments, which are pressurized by a medium with 
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Fig. A.l. Heating rates of the approximation for cores em- 
bedded in spherical filaments, given as function of the ex- 
tinction Ay(0) through the filament center. The assumed 
external pressure is p ex t/k = 2 x 10 4 K/cm 3 . The heating 
rates correspond to the considered ISRFs without (black 
curves) and with dust emission from the ISM (grey lines). 
The heating rates obtained for an ISRF without dust emis- 
sion from the Galactic disc are compared with the accurate 
ray-tracing calculation (symbols). The total heating rates 
(solid lines, filled diamonds) are separated into the sev- 
eral components of the radiation: attenuated external radia- 
tion field, which includes the heating by scattered radiation 
(long dashed lines, open squares); heating by dust emission 
from the filaments (dashed-three dotted curves, open dia- 
monds); heating by dust emission from the cores (dashed- 
single dotted curves, open triangles); and dust emission 
from the Galactic disc (short dashed curve). 



p ext /fc = 2 x 10 4 K/cm 3 



As shown, the simplification for 
leads to highly accurate 
The 



the scattered radiation (Eq. IA.1 
heating rates of the scattered light inside the core 
heating by the dust emission from the filament and the 
core is considerably lower in this approximation because of 
the absence of PAH emission. The difference is more promi- 
nent for the dust heating from the filaments, also because 
the approximation provides only a lower limit. Overall the 
heating by self- absorption inside the cores plays only a mi- 
nor role in the total heating rate. The heating caused by the 
dust emission from the filaments provides a more important 
contribution, which gradually increases with extinction Ay. 
On the basis of the ray-tracing calculation, this contribu- 
tion at Ay = 64 mag is approximately half of the total 
heating rate. However, as the figure suggests, for highly 
embedded cores the dust emission of the galactic disc, in 
particular the emission from PAH molecules, possibly be- 
comes the predominant heating source. 

The effective temperature of the dust emission spectrum 
is obtained in a procedure often used in radiative transfer 
problems where it is assumed that grains at a certain loca- 
tion inside the cloud all have the same temperature. The 
basic difference here is the replacement of the dust emissiv- 
ity by its effective value to obtain the correct dust tempera- 
ture in cores. For the filaments, the intrinsic dust emissivity 
is used. The equations are given here for completeness. The 
effective temperature is determined by the equilibrium of 



the heating and cooling rates 



L hcat = M dust / dA rcf s 4tt j£ eat 



(A.12) 



= M dust J d\Kr^B x {T dust ) = L 

The cooling rate can be written as 
L dust = Mdust 4crgB (« em (^dust)) T^ ust , 



dust 



(A.13) 



where csb is the Stefan-Boltzmann constant. The Planck 
average emission coefficient is given, for a power law (k^ 
^(A/Aq)-' 3 ), by 



„cm 
"A 



/ ™, T u _ «o r(4 + /?)C(4 + /3) fkT dust 

{ ( dustj) " \^ r(4)c(4) 



(A.14) 



where k, h, and c are the Boltzmann constant, the Planck 
constant, and the velocity of light in the vacuum and where 
T(x) and ((x) are the T-function and the ^-function. 



Appendix B: Effective dust properties of the clouds 

In the approximation of the radiative transfer, the absorbed 
flux inside the core or the filament is derived using an ef- 
fective value of the dust extinction and an albedo that is 
based on escape probabilities. To obtain a more accurate es- 
timate of the dust temperature, the values are consistently 
based on the physical density profile of the clouds rather 
than assuming that the dusty medium is homogeneously 
distributed within the given cloud volume. The column 
density related to the corresponding homogeneous sphere 
and cylinder is given in Fig. IB. II In the limit of high over- 
pressure p = Pc/Pcxt, the mean column densities through 
spheres and cylinders are given approximately by 



(Nn) (0) 



Pext/k 



2 x 10 4 K/cm 3 

6.50 x 10 21 p°ciTr 2 
4.33 x 10 21 i5V4 cm- 



sphere), 
2 (cylinder) 



(B.l) 



B.l. Attenuation 

The total attenuation described in terms of an effective 
optical depth caused by the cloud is given by the ratio of the 
flux density entering the cloud from outside by an isotropic 
radiation field to the flux density leaving the cloud without 
being scattered and absorbed. The effective optical depth 
is therefore given by 



-off 



= - In 







-All 







i J dnn F -n,,{n) e" rA(n) | , (B.2) 



where n F is the normal vector to the surface element the 
flux density is derived for, n / is the vector of the photons 
illuminating the surface element, and t\(CI) the total opti- 
cal depth in direction J7. For a self-gravitating cloud, the 
attenuation depends (apart from the actual shape) on the 
density profile determined by the overpressure p = p c /.Pext 
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Fig. B.l. Column density at the center of a self-gravitating 
sphere or a cylinder as a function of the over pressure 
p = Pc/Pcxt, where the mass is homogeneously distributed 
within the cloud radius. The column density and the extinc- 
tion are given for an external pressure p cx t = 2 x 10 4 K / cm 3 . 
The asymptotic behavior is shown as dashed lines. The 
symbols mark the positions of the accurate radiative trans- 
fer calculations. For comparison also the true column den- 
sities at the cloud centers is shown (light grey curves). 



and the radial optical depth t a . For the sphere, the effective 
optical depth is 



71-/2 



T x (p,t\) = -In ^ 2 / ill) sill l) a,:<l>< 



and for the cylinder 



-r A (sini9) 



(B.3) 



tt/2 

(j5,t a ) = -lni 2 J dtf sin?? 
^ o 

tt/2 

x J dip sin cosipe~ T ^ sinlp)/sinm L 
o ' 

where the optical depth in the exponent is given by 



(B.4) 



only the contributions close to the rim, the effective optical 
depth of a sphere becomes 



rf(p,r A »l) a 21n{2i? c i^ x 7^Hn H (-Rci)} 

In 2, sphere 

3 

4 ' 



In I, cylinder. 



(B.7) 



The effective extinction of a sphere and a cylinder for a wide 
range of overpressures are given in Fig. IB.2I For high over- 
pressure p, the effective optical depth can be considerably 
lower than the corresponding mean optical depth through 
a comparable homogeneous cloud of the same mass. 



B.2. Escape probability 

The escape probability is based on the assumption that the 
emissivity is proportional to the density. It is given by the 
ratio of the luminosity escaping the cloud without being 
scattered or absorbed to the total luminosity produced in 
the cloud 



Pesc(P,T A ) 



where 

F-(p, 0) = V cl 4n (W«h) («h> /Sol 
and 

F~(p,r x ) = Jdttn F nj J dh]\(l,tt) 
xexpj-^ d/'A$ xt (/',ft) j . 
For the sphere, the escape probability becomes 
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(B.10) 
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and for the cylinder 
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J dl n H (R c i\/l 2 + sin 2 0) e~ T ^ l ' s ^Xl) 
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tt/2 tt/2 
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Since for self-gravitating clouds, R c i cx T j ^/p cx t and jih oc 
Pcxt/T, where T is the effective temperature of the gas in 
the cloud that combines the different velocity components 
the extinction depends for a given overpressure only on the 
square root of the external pressure (t a oc -Jpext)- 

At low optical depths, the effective optical depth be- 
comes equal to the mean optical depth 

rf{p, r A « 1) = (r A > = i(V cl /S cl ) (n H ) (ftf/fin). (B.6) 

In case of highly opaque clouds, most of the light contribut- 
ing to the flux density arises from the rim where the 
extinction is lowest. The optical depth close to the edge 
is t\ R c \ (fe A xt /W) nn(R c {) 2(n/2 — <p). By considering 



COS V 

x J dln H (R cl ^Jp + sin 2 tp) e-^^^lli) 



— COS $ 

where the optical depth in the exponent is given by 

I 



T\(l,s'mi3) — i? c i(K A xt /riH) / dZ'n H (iiciA// /2 + sin 2 0).(B.13) 

— cos # 

The escape probability for self-gravitating spheres and 
cylinders is shown in Fig. IB.31 For intermediate values, the 
escape probability is lower than for the homogenous clouds. 
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Fig. B.2. Effective optical depth of self-gravitating spheres and cylinders as function of the mean optical depth for 
several assumptions for the over pressure p — p c /Pcxt- The overpressure is varied from 10° to 10 3 in steps of A logp = 0.5. 
The dashed line shows the effective optical depth for clouds with over pressure p = 14.04, the over pressure for critical 
stable spherical clouds. The thin grey lines show the approximations for low and high mean optical depths. The effective 
extinction of the cylinders is also compared to the corresponding values of a homogeneous sphere of the same mean 
optical depth (dotted line). 
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Fig. B.3. Escape probability of self-gravitating spheres and cylinders as a function of the mean optical depth. The 
overpressure is varied in the same range as in Fig. IB.2I The dashed line shows the escape probability of a critical stable 
cloud with p = 14.04 and dotted line in the right-hand figure the escape probability of a homogeneous sphere. 



B.3. Cloud albedo 

The cloud albedo is based on the escape probability de- 
scribed in the previous section. It is assumed that the scat- 
tered light is produced in proportion to the density. For 
multiple scattered events, this gives the cloud albedo 



U\P csc (p,T\) 
1 - U>\(1 -p esc (p,T\)) ' 



(B.14) 



These values are certainly very approximately. However, 
the formula indicates to first order the general way in which 
the cloud albedo decreases for higher extinction values. As 
a consequence, a larger fraction of light is absorbed inside 
the cloud. 



References 

Arendt, R. G., Odegard, N., Weiland, J. L., et al. 1998, ApJ, 508, 74 
Bernard, J. P., Boulanger, F., Desert, F. X., & Puget, J. L. 1992, 
A&A, 263, 258 



Bernard, J. P., Boulanger, F., & Puget, J. L. 1993, A&A, 277, 609 

Bohlin, R. C, Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 

Boisse, P. 1990, A&A, 228, 483 

Bok, B. J. & Reilly, E. F. 1947, ApJ, 105, 255 

Bonnor, W. B. 1956, MNRAS, 116, 351 

Boulares, A. & Cox, D. P. 1990, ApJ, 365, 544 

Chakrabarti, S. & McKee, C. F. 2008, ApJ, 683, 693 

Curry, C. L. & McKee, C. F. 2000, ApJ, 528, 734 

Doty, S. D., Metzler, R. A., & Palotti, M. L. 2005, MNRAS, 362, 737 

Draine, B. T. & Li, A. 2007, ApJ, 657, 810 

Dunne, L. & Eales, S. A. 2001, MNRAS, 327, 697 

Ebert, R. 1955, Zeitschrift fur Astrophysik, 37, 217 

Evans, II, N. J., Rawlings, J. M. C, Shirley, Y. L., & Mundy, L. G. 

2001, ApJ, 557, 193 
Ferriere, K. M. 2001, Reviews of Modern Physics, 73, 1031 
Fiege, J. D. & Pudritz, R. E. 2000a, MNRAS, 311, 85 
Fiege, J. D. & Pudritz, R. E. 2000b, MNRAS, 311, 105 
Fischera, J. 2004, A&A, 428, 99 
Fischera, J. & Dopita, M. 2004, ApJ, 611, 911 
Fischera, J. & Dopita, M. 2005, ApJ, 619, 340 
Fischera, J. & Dopita, M. 2008, ApJS, 176, 164 
Fischera, J., Dopita, M., & Sutherland, R. 2003, ApJL, 599, L21 
Fitzpatrick, E. L. 1999, PASP, 111, 63 



24 



Jorg Fischera: On the SEDs of passively heated condensed 



Gordon, K. D., Misselt, K. A., Witt, A. N., & Clayton, G. C. 2001, 

ApJ, 551, 269 
Hcgmann, M. & Kegel, W. H. 2003, MNRAS, 342, 453 
Kauffmann, J., Pillai, T., Shetty, R., Myers, P. C, & Goodman, A. A. 

2010, ApJ, 716, 433 
Kruegel, E. & Siebenmorgen, R. 1994, A&A, 288, 929 
Larson, R. B. 1981, MNRAS, 194, 809 
Larson, R. B. 1985, MNRAS, 214, 379 

Mathis, J. S., Mezger, P. G., & Panagia, N. 1982, A&A, 105, 372 
Netterfield, C. B., Ade, P. A. R., Bock, J. J., et al. 2009, "ApJ", 707, 
1824 

Nutter, D., Stamatcllos, D., & Ward-Thompson, D. 2009, MNRAS, 
396, 1851 

Olmi, L., Ade, P. A. R., Angles-Alcazar, D., et al. 2009, "ApJ", 707, 
1836 

Ossenkopf, V. & Henning, T. 1994, A&A, 291, 943 
Ostriker, J. 1964, ApJ, 140, 1056 

Stamatellos, D. & Whitworth, A. P. 2003, A&A, 407, 941 
Stamatellos, D., Whitworth, A. P., Andre, P., & Ward-Thompson, D. 

2004, A&A, 420, 1009 
Stodolkicwicz, J. S. 1963, Acta Astronomica, 13, 30 
Varosi, F. & Dwek, E. 1999, ApJ, 523, 265 
Weingartner, J. C & Draine, B. T. 2001, ApJ, 548, 296 
Witt, A. N. & Gordon, K. D. 1996, ApJ, 463, 681 
Witt, A. N. & Gordon, K. D. 2000, ApJ, 528, 799 
Wolf, S., Fischer, O., & Pfau, W. 1998, A&A, 340, 103 



